C.  A.  Blackmon,  N.  R.  Murphy,  Jr. 


TECHNICAL  REPORT  NO.  3-783 


Reproduced  bv 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

Springfield.  Va  221SI 


AN  ANALYTICAL  MODEL.  FOR  PREDICTING 
CROSS-COUNTRY  VEHICLE  PERFORMANCE 

APPENDiX  C:  VEHICLE  PERFORMANCE  IN  VERTICAL 
OBSTACLES  (SURFACE  GEOMETRY) 

b,  D  D  C 


0 


* 


February  1972 


Sponsored  by  Advanced  Research  Projects  Agency  and  Directorate  of  Research, 
Development  end  Engineering,  U.  S.  Army  Materiel  Command 


Samic*  Agency  U.  S.  Army  Materiel  Command 


Conducted  by  U.  S.  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  Mississippi 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


i 


IS*. 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  THE  BEST 
QUALITY  AVAILABLE. 

COPY  FURNISHED  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


*  F 


The  findings  in  this  report  ore  not  to  be  construed  os  on  official 
Department  of  the  Army  position  unless  so  designated 
by  other  authorised  documents. 


Destroy  this  report  when  no  longer  needed.  Do  not  return 
it  to  the  originator. 


BISTRES !!im  CCB63 


BIST.  I  AVAIL  tni'X  SFTCIAl 


TECHNICAL  REPORT  NO.  3-783 


AN  ANALYTICAL  MODEL  FOR  PREDICTING 
CROSS-COUNTRY  VEHICLE  PERFORMANCE 

APPENDIX  C:  VEHICLE  PERFORMANCE  IN  VERTICAL 
OBSTACLES  (SURFACE  GEOMETRY) 

by 

C.  A.  Blackmon,  N.  R.  Murphy,  Jr. 


D  D  C 


TO 


n  n 


MAR  6  TO 

I 


February  1972 

Sponsored  by  Advanced  Research  Projects  Agency  and  Directorate  of  Research, 
Development  and  Engineering,-  U.  S.  Army  Materiel  Command 

Service  Agency  U.  S.  Army  Materiel  Command 
Project  Nos.  I-T-0-62I I2-A-I3I  and  I -T-0-62 1 03- A-046-02 

Conducted  by  U.  S.  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  Mississippi 

ARMY-MRC  VICK9BURQ,  MIBB 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


Unclassified _ 

1 

.  „  ,  ,  DOCUMENT  CONTROL  DATA  -  R  &  0 

J<  W  -W  .H.  . 


v  (Corpormlt  muthoe) 

u.  S.  Army  Engineer  Waterways  Experiment  Station 
Vieksburg,  Mississippi 


""0»T  6CCU.ITV  CL.MIFICATION 

Unelassified 


16.  G.OU. 


AN  ANALYTICAL  MODEL  FOR  PREDICTING  CROSS-COUNTRY  VEHICLE  PERFORMANCE-  »pnPmnT  „ 
VEHICLE  PERFORMANCE  IN  VERTICAL  OBSTACLES  ( SURFACE  SrSf  ’ 


|*  OI.CPIPTIVI  NOTE!  (Typ,  ol  ,^c„  w  Inclu.l,.  dp,..) - - - 

^Appendix  C  to  a  final  report  not  yet  published 


|»-  *u  TMO.Iil  (M,t.  t  Mn>,  mlddlm 


Claude  A .  Blackmon 
Newell  R,  Murphy,  Jr. 


|S>  "IPORT  OATI 

February  1972 


P  CONTRACT  OR  CHANT  NO.  - " - 

ARPA  Order  Ho.  400 

*.  P  NOJEC  T  NO. 

*•  1-T-0-62112-A-131  and 

1-T -0-62103 -A  -01(6-02 

|To.  OI.TRI  BUT  (ON  tMTIWHT - 

Approved  for  public  release;  distribution  unlimited. 

• lU^MLIMINTART NOTE!  - - - 


|7*.  TOTAL  NO.  OF  PAGES 

ZjlO— 


|7fc.  NO.  or  WE  PI 

26 


**•  ORIQINATON'9  NIRORT  NUfcfBE  H(9)  1 

Technical  Report  No.  3 -783)  Appendix  C 


*■  S/T.Mij,r,"°"T  no'*’  <*nr  6. 


II.  ..ONIONING  MILITARY  ACTIVITY  - - - - 

Advanced  Research  Projects  Agency  and 
Directorate  of  Research,  Development  and 
Engineering,  U.  g.  Army  Materiel  Command 


and  MR7  t™,v  1  dynamic  response  prediction  models  for  the  M60A1  tank 

peak  vertical '  ’’  a"d  ^re'Jicted  vehicle  performance  in  terms  of 

peaK  vertical  t  1  longitudinal  accelerations  for  78  rigid  ob-'ar-ie  l+vl  i 

-  •»“ « »  ssmss  zirz  rfe  r;tr  refi— 


DD  i  n*v**.1473  »■*»*■*’ SoiTSSy 


I  JAM  94,  WHICH  19 


Unclassified 

fpcurity  C'la  tiiricallon 


Unclassified 

Security  CliMtricutton 


THE  CONTENTS  OF  THIS  REPORT  ARE  NOT  TO  BE 
USED  FOR  ADVERTISING,  PUBLICATION,  OR 
PROMOTIONAL  PURPOSES.  CITATION  OF  TRADE 
NAMES  DOES  NOT  CONSTITUTE  AN  OFFICIAL  EN¬ 
DORSEMENT  OR  APPROVAL  OF  THE  USE  OF  SUCH 
COMMERCIAL  PRODUCTS. 


FOREWORD 


The  study  reported  herein  was  performed  by  the  U.  S.  Army  Engi¬ 
neer  Waterways  Experiment  Station  (WES)  for  the  Office,  Secretary  of 
Defense  (OSD),  Advanced  Research  Projects  Agency  (ARPA),  and  is  a  por¬ 
tion  of  one  task  of  the  overall  Mobility  Environmental  Research  Study 
(MERS)  sponsored  by  OSD/ARPA  for  which  the  WES  was  the  prime  contractor 
and  the  U.  S.  Army  Materiel  Command  (AMC)  was  the  service  agent.  The 
broad  mission  of  Project  MERS  was  to  determine  the  effects  of  the  var¬ 
ious  features  of  the  physical  environment  on  the  performance  of  cross¬ 
country  ground  contact  vehicles  and  to  provide  therefrom  data  that  can 
be  used  to  improve  both  the  design  and  employment  of  such  vehicles.  A 
condition  of  the  project  was  that  the  data  be  interpretable  in  terms  of 
vehicle  requirements  for  Southeast  Asia.  The  funds  employed  for  this 
study  wore  allocated  to  WES  through  AMC  under  ARPA  Order  No.  400.  Some 
funds  for  preparation  and  publication  of  this  report  were  provided  by 
the  Research,  Development  and  Engineering  Directorate,  AMC,  under  Depart¬ 
ment  of  the  Army  Project  1-T-O -62112 -A-131,  "Environmental  Constraints 
on  Materiel,"  and  Project  1-T-0-62103-A-046-02,  "Surface  Mobility." 

This  appendix  is  one  of  seven  to  a  report  entitled  An  Analytical 
Model  for  Predicting  Cross-Country  Vehicle  Performance  (in  preparation). 
These  appendixes  are: 

A.  Instrumentation  of  Test  Vehicles  (July  1967) 

B.  Vehicle  Performance  in  Lateral  and  Longitudinal  Obstacles 
(Vegetation) 

Volume  I:  Lateral  Obstacles  (December  1968) 

Volume  II:  Longitudinal  Obstacles  (July  1968) 


v 


C.  Vehicle  Performance  in  Vertical  Obstacles  (Surface 
Geometry)  (February  1972) 

D.  Performance  of  Amphibious  Vehicles  in  the  Water-Land 
Interface  (Hydrologic  Geometry)  (February  1970) 

E.  Quantification  of  the  Screening  Effects  of  Vegetation  on 
Driver's  Vision  and  Vehicle  Speed  (April  1971) 

F.  Soil-Vehicle  Relations  on  Soft  Clay  Soils  (Surface  Compo¬ 
sition)  (August  1970) 

G.  Application  of  Analytical  Model  to  United  States  and 
Thailand  Terrains  (in  preparation) 

This  study  was  conducted  by  personnel  of  the  Mobility  and  En¬ 
vironmental  (M&E)  Division,  under  the  general  supervision  of  Mr.  W.  J. 
Turnbull,  Technical  Assistant  for  Soils  and  Engineering  (retired); 

Mr.  W.  G.  Shockley  and  Mr.  S.  J.  Knight,  Chief  and  Assistant  Chief,  re¬ 
spectively,  M&E  Division;  and  under  the  direct  supervision  of  Mr.  A.  A. 
Hula,  Chief,  Vehicle  Studies  Branch  (VSB). 

The  field  tests  reported  herein  were  conducted  during  the  period 
January -June  1968  by  Mr.  J.  L.  Gargaro,  VSB,  and  G.  Switzer,  Mobility 
Research  Branch  (MRB);  the  mathematical  models  of  the  M60A1  tank  and 
M37  truck  used  to  make  the  performance  predictions  were  formulated  and 
described  by  Mr.  N.  R.  Murphy,  Jr.,  MRB,  and  the  analysis  of  the  data 
presented  was  performed  by  Mr.  C.  A.  Blackmon,  VSB.  This  report  was 
prepared  by  Messrs.  Blackmon  and  Murphy. 

Acknowledgment  is  made  to  Mr.  W.  A.  Gross,  Jr.,  and  Mr.  P.  R. 

Gula,  Development  and  Proof  Services,  Aberdeen  Proving  Ground  (APG), 
for  their  support  and  assistance  during  the  conduct  of  the  field  tank 
tests  reported  herein. 

Directors  of  the  WES  during  this  study  and  the  preparation  of  this 
report  were  COL  Alex  G.  Sutton,  Jr.,  CE,  COL  John  R.  Oswalt,  Jr.,  CE, 

COL  Levi  A.  Brown,  CE,  and  COL  Ernest  D.  Peixotto,  CE.  Technical  Direc¬ 
tors  were  Messrs.  J.  B.  Tiffany  and  F.  R.  Brown. 
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CONVERSION  FACTORS,  BRITISH  TO 


British  units  of  measurement  used  in 
metric  units  as  follows: 


square  inches 
cubic  feet 


pounds 

kips 

short  tons  (2000  lb) 
pounds  per  square  inch 
foot-pounds 
miles  per  hour 


2. 5i+ 
0.301+8 
6.1+516 
0.0283168 
O.U5359237 
453.5923  7 
907.185 
0.070307 
0.138255 


1.609344 
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UNITS  OF  MEASUREMENT 


is  report  can  be  converted  to 

_ To  Obtain _ 

centimeters 

meters 

square  centimeters 
cubic  meters 
kilograms 
kilograms 
kilograms 

kilograms  per  square  centimeter 
meter -kilograms 
kilometers  per  hour 


SUMMARY 


This  appendix  presents  a  brief  history  of  vehicle  dynamics  mod- 
eling,  a  recapitulation  of  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  approach  to  the  problem  of  predicting  vehicle  performance  in 
terrain  containing  discrete  vertical  obstacles,  and  descriptions  of  dy¬ 
namic  response  prediction  models  for  the  M60A1  tank  and  K37  truck  and 
compares  measured  and  predicted  vehicle  performance  in  terms  of  peak 
and  longitudinal  accelerations  for  78  rigid  obstacle  tests 
with  the  two  vehicles.  Major  conclusions  from  the  tests  were  that  per¬ 
formances  of  the  M60A1  tank  and  M37  truck  in  terms  of  peak  vertical  and 
peak  longitudinal  accelerations  experienced  at  the  driver's  seat  when 
raversing  discrete,  rigid  obstacles  can  be  correlated  with  impact 
speed,  that  the  mathematical  techniques  described  yield  reasonably  ac- 
cura^e  predictions  of  the  speed  at  which  the  M60A1  tank  can  contact  a 
rigid  obstacle  without  exceeding  specified  tolerance  limits,  and  that 

refinement  is  needed  in  the  dynamic  response  prediction  model  for  the 
Mo  (  truck. 
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AN  ANALYTICAL  MODEL  FOR  PREDICTING 


CROSS-COUNTRY  VEHICLE  PERFORMANCE 


APPENDIX  C:  VEHICLE  PERFORMANCE  IN  VERTICAL 
OBSTACLES  (SURFACE  GEOMETRY) 


PART  I:  INTRODUCTION 


Background 

l •  The  main  text  of  this  report  (not  yet  published)  describes  the 
development  of  an  analytical  model  for  predicting  the  cross-country  per¬ 
formance  of  a  vehicle.  The  model  was  based  on  an  energy  concept  within 
the  framework  of  classical  mechanics  which  requires  that  cause -and - 
effect  relations  be  established  between  discrete  terrain  factors  and  ve¬ 
hicle  response.  This  appendix  deals  with  the  effects  of  a  single  ter¬ 
rain  factor- -vertical  obstacles.  The  term  "obstacle"  in  general  refers 
to  all  features  of  the  terrain,  except  soil,  that  are  inhibitory  to  ve¬ 
hicle  mobility.  The  obstacle -effects  spectrum  of  vehicle  mobility 
ranges  from  complete  immobilization  to  minor  speed  reduction.  For  the 
purpose  of  the  overall  study,  obstacles  were  categorized  according  to 
the  direction  of  motion  forced  upon  a  vehicle  negotiating  the  obstacle, 

i.e.  vertical,  lateral,  or  longitudinal. 

2.  Perhaps  the  most  universal  single  terrain  feature  that  pro¬ 
duces  an  inhibiting  effect  on  vehicle  ground  mobility  is  small-scale 
surface  geometry.  Surface  geometry  features  occur  in  a  bewildering  ar¬ 
ray  of  sizes  and  configurations,  and  produce  effects  on  vehicles  that 
range  from  "vibration"  to  "shock"  to  "immobilization,"  depending  on  the 
speed  and  size  of  the  vehicle  in  relation  to  the  size  and  spacing  of 
obstacles. 

3.  Vibration-producing  features  are  those  surface  irregularities 
of  heights  that  can  be  measured  in  inches*  rather  than  feet,  and  of 

*  A  table  of  factors  for  converting  British  units  of  measurement  to 
metric  units  is  presented  on  page  ix. 
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distance  between  features  that  can  be  measured  in  feet  rather  than 
tens  of  feet.  The  dominant  interaction  for  vibration -producing  fea¬ 
tures  is  dynamic  and  is  concentrated  in  the  action  of  the  vehicle  sus¬ 
pension  as  the  vehicle  passes  over  the  features  at  a  speed  in  excess  of 
creep  speed,  say  5  mph  or  more,  with  the  vehicle  speed  limited  by  con¬ 
siderations  of  driver  and/or  cargo  safety. 

4.  Shock -inducing  features  occur  as  discrete  "bumps"  and  produce 
a  dynamic  interaction  that  also  is  concentrated  in  the  vehicle  suspen¬ 
sion  system.  The  dynamic  action  associated  with  these  features  is  of  a 
high-amplitude,  low-frequency  nature.  The  vehicle  speed  is  ultimately 
linuted  by  considerations  of  driver  and/or  cargo  safety. 

5.  Features  that  are  likely  to  produce  i:nmobilizations  occur  as 
discrete  obstacles,  as  do  the  shock -inducing  features;  however,  the 
former  are  of  such  size  and  shape  that  their  negotiation  can  be  consid¬ 
ered  in  terms  of  static  phenomena  since  an  attempt  to  traverse  them 
must  be  made  at  creep  speed.  The  controlling  factors  of  obstacle - 
vehicle  interaction  for  these  features  are  the  geometry  of  the  feature, 
the  geometry  of  the  vehicle  body  and  running  gear,  and  the  ability  of 
the  vehicle  to  exert  sufficient  tractive  effort  to  lift  itself  over 
the  feature. 

6.  Early  testing  did  little  more  than  outline  the  problem.  It 
was  apparent  that  nothing  less  than  an  elaborate  computer  program  would 
suffice,  hence  the  major  effort  of  this  study  was  directed  toward  exam¬ 
ination  of  current  mathematical  modeling  techniques  and  the  refinement 
thereof. 


Purpose  and  Scope 

7.  This  appendix  presents  a  brief  history  of  vehicle  dynamics 
modeling,  a  discussion  of  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  approach  to  the  problem  of  predicting  vehicle  performance 
when  traversing  discrete  vertical  obstacles,  a  description  of  the  dy¬ 
namics  submodels  used,  and  a  comparison  of  measured  and  predicted  ve¬ 
hicle  performance . 
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8.  Seventy-eight  obstacle  -vehicle  tests  were  run  with  two  ve¬ 
hicles  at  several  speeds  over  a  range  of  obstacle  heights.  Dynamic  re¬ 
sponse  predictions  in  terns  of  peak  vertical  acceleration  and  peak  lon¬ 
gitudinal  acceleration  were  made  for  34  combinations  of  obstacle  height 
and  vehicle  speed.  The  tests  and  predictions  were  limited  to  the  ve¬ 
hicles  crossing  a  rigid,  nondeformable  obstacle  on  a  uniformly  hard 

surface . 


A  History  of  Vehicle  Dynamics  Modeling 


In  retrospect 

9.  Since  the  advent  of  the  automobile  (especially  since  the 
1920's),  research  by  government  and  industry  has  been  conducted  on  a 
continuing  basis  in  various  countries  with  regard  to  highway  design  and 
construction  and  also  vehicle  design,  with  particular  emphasis  on  steer¬ 
ing  control,  power  train,  and  suspension-the  principal  contributors  to 
the  safety,  efficiency,  and  riding  comfort  of  on-road  vehicles. 

10.  The  speed  at  which  a  driver  of  a  vehicle  will  traverse  ob¬ 
stacles  or  continuous  irregular  terrain  is  controlled  primarily  by  the 
level  of  vibration  activity  that  does  not  exceed  his  particular  ride 
comfort  level.  Vehicle  vibration  or  ride  is  sensed  by  a  driver  or  pas¬ 
senger  through  sight,  touch,  and  hearing  in  response  to  external  stim¬ 
uli,  such  as  motions,  forces,  and  sounds.  Whenever  this  sensation  be¬ 
comes  too  severe,  the  driver  will  alter  the  vehicle's  speed  until  the 
sensation  reaches  an  acceptable  level.  This  sensation,  therefore,  is  a 
significant  factor  in  determining  the  speed  of  a  vehicle  over  a  gi.en 
terrain.  The  irregular  terrain -vehicle  problem  is  essentially  one  of 
dynamics,  and  its  solution  must  include  the  combined  effects  of  the 
surface  being  traversed,  the  vehicle,  and  the  driver.  Because  of  the 
complexity  of  the  problem  and  the  desire  to  produce  better  riding  ve¬ 
hicles,  considerable  effort  has  been  expended  on  modeling  dynamic  ve- 

hide  response.  . 

11.  Because  of  the  lack  of  mathematical  techniques  required  in 

modeling  suspension  systems,  much  of  the  early  work  consisted  of 
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cut -and -try  methods.  The  first  significant  contributions  to  an  analyt¬ 
ical  treatment  of  vehicle  dynamics  were  performed  by  Rowell,1  Guest,2 
and  Olley3  in  the  early  1920's  and  30 's. 

12.  In  1941,  a  mechanical  differential  analyzer  was  built  by 
Schilling  and  Fuchs**  specifically  for  suspension  analysis,  and  although 
it  was  suited  to  only  a  single-degree -of -freedom  system,  it  did  permit 
the  inclusion  of  the  nonlinear  characteristic  of  shock  absorbers.  The 
analyzer  was  used  in  the  continuous  determination  of  transient  motions 
and  in  the  portrayal  of  the  effect  on  motion  by  changes  in  the  charac¬ 
teristics  of  the  shock  absorber.  This  differential  analyzer  was  the 
forerunner  of  today's  analog  computer,  and  its  development  led  to  rapid 
advances  in  suspension  analysis  and  design.  By  the  1950's,  it  was 
widely  exploited  by  the  automotive  industry. 

13-  In  1953,  Jeska  developed  a  four -degree -of -freedom  model 
that  included  pitch  and  bounce  of  the  body  and  vertical  motions  of  the 
front  and  rear  wheels.  The  forcing  function  was  an  actual  road  wave 
measured  by  a  photographic  technique.  In  1955,  Bodeau,  Bollinger,  and 
Lipkin  of  Ford  Motor  Company  developed  a  detailed  ride  analysis  in 
which  a  nine -degree -of -freedom  model  was  used  to  describe  a  passenger 
car.  In  i960,  Kohr^  of  General  Motors  Corporation  developed  a  mathe¬ 
matical  simulation  of  automobile  ride.  In  his  simulation,  a  measured 
road  profile  was  recorded  on  magnetic  tape,  and  the  tape  was  fed 
through  an  analog  computer  model  of  the  vehicle  to  predict  the  vehicle 
motions,  i.e.  pitch,  bounce,  and  roll.  The  resulting  motions  were  used 
to  drive  a  vibration  simulator,  which  was  used  as  a  laboratory  means  of 
assessing  the  (-ffect  of  the  vibration  on  humans. 

14.  Until  about  i960,  the  analysis  of  ride  had  been  concerned 
primarily  with  the  suspension  s^  stem  and  means  of  improving  the  ride 
quality.  Although  considerable  work  was  done  in  the  area  of  human  tol¬ 
erance  to  vibration,  a  means  for  quantifying  human  tolerance  to  vibra¬ 
tions  had  not  been  developed.  Van  Deusen^  has  shown  that  very  little 
of  the  research  done  actually  pertains  to  the  off-road  environment. 

Most  experiments  have  been  devised  to  assess  human  response  to  sinusoi¬ 
dal  motion  in  only  one  direction,  while  the  more  complex  ride  comfort 
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problems  involve  random  vibrations  in  various  directions.  The  most 
frequently  used  criteria  have  been  those  of  Diectanan9  and  Janeway,10 
who  developed  simple  formulas  for  relating  comfort  limits  to  amplitude 
and  frequency  of  vibration.  There  have  been  several  studies  of  "on- 
the-road"  measurements  of  ride  comfort.  For  example,  Von  Eldik  Thieme11 
examined  the  Dieckman -Janeway  criteria  in  the  actual  vehicle  environ¬ 
ment,  but  he  met  with  little  success.  Van  Deusen18’13  and  Versace14 
used  a  technique,  referred  to  as  cross  modality,  in  which  subjects  re¬ 
ceived  noise  signals  through  earphones  and  adjusted  and  matched  the 
signal  s  level  to  the  sensation  level  of  ride  vibration,  a  statistical 
analysis  showed  favorable  correlations  of  the  measured  accelerations 
With  ride  sensation,  and  at  least  indicated  that  correlations  between 
ride  sensation  and  vibration  were  possible. 

15.  In  the  late  1950's  the  Department  of  Defense  began  to  recog¬ 
nise  the  significance  of  vehicle  vibration  on  off-road  mobility  and 
-pen ^efficiency.  An  extensive  effort  was  begun  to  quantify  the  vehi- 
=le  VI  ration  problem  and  to  chelate  it  with  hunmn  response  and  ter- 
rain  characteristics.  Interest  was  shifts  r  . 

chastic  techniques  The  3  T  “  -to- 

■  ter  teehr;la-'Je  consists  of  classifying  terrain 

tistic  n  ”1 '  Pertlnel,t  StatlStICS  M1  analy2l"g  sta- 

stically.  The  groundwork  for  this  type  of  analysis  was  begun  in  1959 

by  Bogdanoff  and  Kozin,19  who  described  in  detail  the  statist^!  anal! 

ysis  of  the  responses  of  simple  linear  systems  to  random  terrain  inputs. 

Although  the  vehicle  mode!  s  were  simple  and  idealised,  the  analyses 

provided  a  starting  point  and  yielded  much  useful  information  regarding 

fundamental  relations  between  pertinent  vehicle  parameters  and  statis- 

Ilrm  rln  qua"l1ties-  studies  of  Bieniek18 

I960),  Van  Deusen  (1962).  and  Bussman18  (!<«,),  „ho  foUowa  eSEen. 

rally  the  same  approach  as  that  described  by  Bogdanoff  and  Kotin.  The 
one  notable  exception  was  Van  Deusen's  introduction  of  a  nonlinear  ve- 
hide  system  into  his  statistical  analysis. 

1  i.  Organized  discrete  obstacle-vehicle  research  in  the  Western 
world  was  given  special  attention  as  the  result  of  World  War  II  experi¬ 
ences.  Early  U.  S.  military  efforts  were  concerned  with  designing 


military  vehicles  that  would  reduce  immobilizations  caused  by  obstacle 
interference.  Obstacle  test  courses  were  constructed  at  Aberdeen  Prov¬ 
ing  Ground  ( APG) ,  and  tests  on  these  courses  have  become  part  of  the 
overall  vehicle  engineering  evaluation  test  program.  Results  of  these 
studies  led  to  the  recent  development  of  articulated  vehicles.  Dis¬ 
crete  obstacle -vehicle  research  studies  in  the  United  States  gained 
more  emphasis  about  the  mid-1950's  when  terrain  factors  other  than 
soils  were  introduced  as  deterrents  to  off -road  vehicle  travel,  and 
more  attention  was  given  to  obstacle  geometry  interference  and  the  ef¬ 
fects  of  dynamic  response  on  vehicle  performance.  By  the  early  1960's 
these  studies  produced  several  static  and  quasi-dynamic  models  which 
related,  by  simple  two-dimensional  static  mechanics,  slope  and  obstacle 
geometry  to  go-no  go  performance. 

17.  In  1963,  the  U.  S.  Army  Tank -Automotive  Command  (TACOM)  be¬ 
gan  basic  research  on  the  effects  of  vehicle  vibration  on  human  re¬ 
sponse.  This  work  was  based  on  the  results  of  past  studies  and  so  was 
oriented  toward  quantifying  the  effects  of  random  vibration  on  vehicle- 
driver  performance.  Two  performance  parameters  were  developed  to  de¬ 
scribe  human  response — acceleration  density  and  absorbed  power.  Of  the 
two  performance  parameters,  absorbed  power  is  preferred  in  quantifying 
human  response  to  vibration  since  it  is  a  descriptor  of  the  flow  of  en¬ 
ergy  from  the  vibrating  vehicle  to  the  driver.  This  led  to  the  TACOM 
V-ride  concept  in  which  ride  limiting  speed  is  determined  as  that  speed 
at  which  the  driver's  absorbed  power  reaches  6  watts.  During  the 
1960's  the  Department  of  Defense  sponsored  several  studies  in  the  de¬ 
velopment  and  application  of  ad  hoc  comprehensive  cross-country  models 
in  which  V-ride  was  incorporated  as  a  submodel. 

18.  In  the  early  1960's  the  scope  of  WES  mobility  research  was 

expanded,  and  static  and  dynamic  surface  configuration,  vehicle,  driver 

interaction  studies  were  initiated.  In  1965,  FMC  Corporation  conducted 
19 

a  study  for  WES  to  determine  the  feasibility  of  using  a  digital  com¬ 
puter  to  simulate  the  dynamic  response  of  ground  vehicles  traveling 
over  unyielding  irregular  terrain  segments.  This  study  resulted  in  the 
development  of  a  generalized  mathematical  model  of  an  n-axle  vehicle 
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which,  within  limits,  is  suitable  for  both  wheeled  and  tracked  vehicles. 
Current  approaches 


19.  Today's  practice  in  modeling  the  effects  of  surface  config¬ 
uration  on  vehicle  performance  consists  essentially  of  two  types  of 
analysis  which  are  separated  on  the  basis  of  the  kind  of  vehicle  and/or 
driver  interaction  anticipated.  Regardless  of  the  analysis  performed, 
the  terrain  profile  and  associated  discrete  obstacles  generally  are 
considered  rigid.  This  consideration  represents  the  worst  conditions 
that  might  be  encountered  from  the  standpoint  of  the  vehicle's  vibra¬ 
tional  behavior.  If  a  terrain  unit  contains  an  irregular  surface  that 
can  be  easily  overridden  by  a  vehicle  without  inducing  frequent  shock, 
vehicle  performance  is  predicted  by  a  dynamic  model.  In  this  case,  the 
problem  is  commonly  identified  as  surface  roughness,  and  the  profile 
used  is  a  statistically  uniform  surface  profile.  If  a  terrain  unit 
contains  discrete  obstacles  larger  than  those  included  in  the  rough 
terrain  analysis  and  which  are  likely  to  produce  immobilizations,  it  is 
assumed  that  traversing  or  circumventing  the  obstacles  will  be  accom¬ 
plished  at  a  creep  speed  and  the  interactions  are  treated  as  static 
phenomena.  The  models  used  in  such  terrain  situations  are  thus  static 
models.  The  controlling  factors  in  the  relation  of  a  static  model  are 
the  geometry  of  the  obstacle  and  vehicle  configuration.  If  the  dis¬ 
crete  obstacles  that  the  vehicle  must  pass  over  occur  at  wide  spacings 
(e.g.  dikes),  and  they  can  be  overridden  at  speeds  greater  than  creep 
speeds,  a  dynamic  model  is  used  to  determine  the  maximum  override  speed. 

20.  Current  vehicle  dynamic  models  simulate  mathematically  the 
dynamic  response  of  selected  points  within  a  vehicle  (usually  at  the 
driver's  seat  or  in  the  cargo  compartment)  as  it  traverses  discrete  ob¬ 
stacles  or  rough  terrain.  Performance  is  generally  expressed  in  terms 
of  relations  between  speed  and  such  response  quantities  as  absorbed 
power,  root -mean -square  acceleration,  or  peak  acceleration,  and  then 
referenced  to  established  tolerances  for  horizontal  and  vertical  accel¬ 
erations  and  power  limits.  Most  dynamic  models  predict  only  vertical 
motions;  however,  they  can  readily  be  modified  to  include  horizontal 
and,  if  necessary,  lateral  motions.  The  mathematical  techniques 
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involved  in  the  formulation  of  dynamic  models  are  common  to  all  models, 
but  differences  occur  in  the  details  of  representing  the  terrain,  ve¬ 
hicle,  and  driver  limits,  and  in  the  size  of  the  computer  required. 

21.  Mathematical  descriptions  of  vehicle  behavior  in  surmounting 
obstacles  are  a  combination  of  dynamic  and  static  models.  Discrete  ob¬ 
stacles  such  as  rocks,  boulders,  mounds,  scarps,  ditches,  etc.,  are 
usually  first  examined  to  determine  whether  or  noo  there  will  be  spa¬ 
tial  interference  between  the  obstacles  and  the  nonpropelling  vehicle 
structure.  This  examination  may  proceed  in  either  two  or  three  dimen - 
sions:  with  or  without  compliance  of  the  vehicle  running  gear,  suspen¬ 
sion,  or  structure;  with  or  without  compliance  of  the  obstacle  itself;  or 
by  spatial  matching  of  the  vehicle  and  the  obstacle,  either  of  which 
may  be  described  more  or  less  completely.  In  some  instances,  the  vehi¬ 
cle  underside  is  modeled  to  scale,  usually  in  two  dimensions.  In 
others,  the  vehicle  shape  is  idealized  to  the  quantitative  description 
and  location  of  salient  features,  and  matching  is  done  through  the  ap¬ 
plication  of  complex  but  ordinary  trigonometry  and  geometry.  The  latter 
procedure  is,  of  course,  more  suitable  for  computer  use.  Where  the 
number  of  obstacle  configurations  assumed  in  an  area  is  relatively 
small,  however,  the  scale -model  experiments  can  be  conducted  once  and 
for  all,  and  the  results  stored  for  subsequent  computer  consultation  as 
needed.  In  addition  to  static  examination  regarding  obstacle-vehicle 
geometry  interference  aspects,  it  is  necessary  also  to  describe  the  dy¬ 
namic  response  of  the  vehicle  when  traversing  a  given  obstacle  at  a 
given  speed.  It  is  this  problem  upon  which  this  report  is  focused. 
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PART  II:  THE  WES  APPROACH 


Obstacle -Vehicle  Interaction  Categories 

22.  Early  efforts  at  WES  were  directed  toward  developing  instru- 
20 

mentation  and  test  procedures  to  measure  and  record  vehicle  response. 
It  was  apparent  that  a  model  for  predicting  the  effects  of  vertical  ob¬ 
stacles  must  consider  at  least  three  broad  categories  of  obstacle - 
vehicle  interactions. 

Obstacle  geometry  interference 

23.  This  model  was  designed  to  answer  such  questions  as:  Is  it 
possible  for  the  vehicle  to  cross  the  obstacle  without  hanging  up?  Is 
the  obstacle  of  such  size  and  configuration  that  the  vehicle  might  be 
in  danger  of  up-ending?  Is  there  sufficient  traction?  Initially,  a 
vwo -dimensional  scale  model  of  the  vehicle  and  of  the  obstacle  was  used 
to  answer  the  first  two  questions  and  to  determine  a  maximum  attitude 
angle  which  was  compared  with  the  maximum  negotiable  slope  to  determine 
a  simple  go-no  go.  A  procedure  for  determining  maximum  slope  negoti- 
able  is  given  in  Appendix  D  of  this  report.  Subsequently,  a  computer 
program  was  written  that  mathematically  determined  the  same  performance 
parame  ters . 

Speed  limited  by  maneuvering 

24.  A  method  for  determining  whether  or  not  the  vehicle  can 
circumvent  the  obstacle  on  the  basis  of  a  theoretical  parameter,  "area 
denied,"  and  empirical  relations  of  area  denied  and  speed  made  good 
were  developed  to  answer  such  questions  as:  Is  it  possible  for  the 
vehicle  to  circumvent  the  obstacles?  What  reduction  in  speed  is 
brought  about  if  the  vehicle  does  circumvent  the  obstacles?  This 
method  is  described  in  Appendix  B,  Volume  I,  of  this  report.  A  com¬ 
puter  program  was  written  that  utilized  a  somewhat  refined  procedure 
for  computing  area  denied  and  the  effect  of  area  denied  on  vehicle 
speed. 

Speed  limited  by  dynamic  response 

25.  If  the  vehicle  can  traverse  the  obstacle,  the  maximum  safe 


speed  at  which  it  can  cross  the  obstacle  must  be  determined.  To  do 
this,  the  performance  parameters  and  conditions  that  limit  speed  must 
be  known.  A  system  for  predicting  vehicle  dynamic  response  was  being 
developed  concurrently  with  the  conduct  of  the  MERS  vehicle  field  test 
programs,  but  the  input  and  output  requirements  of  the  prediction  sys¬ 
tem  were  not  known  at  the  time  the  test  programs  were  conducted.  Con¬ 
sequently,  the  response  data  taken  during  the  MERS  test  programs  did 
not  entirely  satisfy  the  requirements  of  the  prediction  system  to  the 
extent  that  they  could  be  used  for  comparison  with  predicted  values. 

However,  the  test  results  did  tend  to  confirm  previous  cross-country 
23 

studies,  which  show  that  the  maximum  peak  vertical  acceleration  tol¬ 
erated  by  the  driver  while  trying  to  maintain  a  maximum  safe  speed  was 
approximately  2.5  g's.  An  additional  tolerance  limit  of  2.0-g  peak 

longitudinal  acceleration  was  established  on  the  basis  of  competent  ex- 

24 

perience  and  judgment.  Other  tests  confirmed  that  peak  vertical,  ac¬ 
celeration  and  peak  longitudinal  acceleration  could  be  adequately  meas¬ 
ured  in  vehicle  tests.  Computer  models  for  mathematically  simulating 
vehicle  dynamic  response  were  developed  to  predict  peak  vertical  accel¬ 
eration  and  peak  longitudinal  acceleration  for  rigid-frame  wheeled  and 
tracked  vehicles  crossing  discrete  obstacles.  Descriptions  of  the  mod¬ 
els  are  discussed  in  subsequent  paragraphs. 

Dynamic  Response  Prediction  Model  for  M60A1  Tank 

2 6.  A  mathematical  model  of  the  M60A1  tank  was  developed  to  sim¬ 
ulate  the  dynamic  response  of  the  tank  while  crossing  rigid  obstacles 
perpendicular  to  the  path  of  travel.  The  intent  of  this  model  was  to 
portray  as  nearly  as  possible  the  significant  features  of  the  motions 
of  interest,  namely,  those  in  the  vicinity  of  the  driver,  and  to  pre¬ 
dict  obstacle  limiting  speeds  based  on  some  preselected  tolerance  cri¬ 
terion  at  a  minimum  cost.  The  model  thus  does  not  contain  the  detail 
that  might  be  required,  for  example,  in  designing  a  vehicle  suspension 
system. 

27.  A  schematic  diagram  of  the  system  that  was  modeled  is  shown 
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Fig.  Cl.  Schematic  of  M60A1  tank 


in  fig.  Cl.  This  model  consists  of  nine  degrees  of  freedom  that  include 
the  bounce,  pitch,  and  surge  of  the  center  of  gravity  of  the  main  frame 
and  the  vertical  motions  of  each  of  the  six  bogie  wheels.  Motion  in 
this  context  includes  displacement,  velocity,  and  acceleration.  In  ad¬ 
dition,  the  motions  in  the  vicinity  of  the  driver  are  computed.*  The 
geometry  effects  of  the  bogies  are  represented  by  radially  projecting 
stiff  springs  and  the  track  compliance  by  interconnecting  springs  be¬ 
tween  the  bogies  and  three  "feelers"  appropriately  positioned  in  front 
of  the  first  bogie  and  connected  to  it  by  a  spring. 

28.  The  longitudinal  motion  (X  in  reference  axis,  fig.  Cl)  is 
accounted  for  only  in  the  acceleration  determined  from  the  horizontal 
forces  resulting  from  deflections  of  the  bogie  spring  segments.  The 
horizontal  components  of  the  segment  forces  are  summed  for  each  bogie, 
and  the  horizontal  acceleration  is  obtained  by  dividing  this  summation 
by  the  mass  of  the  tank.  This  method  of  accounting  for  horizontal 


Because  of  certain  geometric  interferences,  it  is  not  possible  to 
jn  accelerometer  at  the  driver  position.  An  accelerometer 
intended  to  measure  the  accelerations  in  the  vicinity  of  the  driver 

£a!\b®  £°nvenifntly  located  at  a  position  1  ft  to  the  right  and  2  ft 
behind  the  center  of  the  driver's  seat. 
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accelerations  is  somewhat  analogous  to  towing  a  vehicle  across  an  ob¬ 
stacle,  always  maintaining  a  constant  velocity,  and  determining  the  in¬ 
creased  towing  force  required  to  tow  the  vehicle  over  the  obstacle  at 
the  given  velocity. 

Equations  of  motion 

29.  The  differential  equations  describing  the  motion  of  this 
system  were  developed  by  first  establishing  an  appropriate  set  of  coor¬ 
dinates  and  sign  convention  and  then  placing  the  system  in  a  displaced 
configuration  such  that  each  coordinate  was  affected.  The  relative 
displacements  of  the  masses  produce  forces  on  each  mass  as  shown  by  the 
free  body  diagram  in  fig.  C2.  The  vehicle  characteristics  used  in 
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predicting  vehicle  dynamic  response  ai'e  given  in  table  Cl.  Using 
Newton’s  second  law  of  motion  and  summing  the  forces  on  the  main  frame 
and  on  each  bogie  led  to  the  series  of  equations  listed  below. 

Forces  on  the  main  frame  (sprung  mass): 

r  6  6  .  . 

Mz  =  -I  2)  +  £  cCAj^Aj^  +  Mg 

r  3  3  .  . 

10  =  -  £  k(z\i)Aiii  cos  e  +  l  ctAj^A^  cos  6 
6  5 

-  X  k(A.  )a‘  i  cos  e  -  V  c  (a,  )Aj  cos  0  +  Q 

i=3  111  ^-3  ill 

6 

Mx  =  T  H. 
i=l  1 

Vertical  forces  on  the  bogies  (unsprung  mass): 

•  • 

Mi'z1  =  k(AL)AL  +  c(Al)Al  +  u1B1  -  u0&0  -  M^g  + 

•  m 

MgZg  =  k(Z^)Z^  +  c(Z^)Z^  -  u1&1  +  u252  -  Mgg  + 

M3z3  =  +  c(^)^  -  n2&2  +  ^383  "  +  v3 

Vu  =  k(Au)Ai  +  c(Al)A4  -  ^353  +  4454  -  V  +  V4 

•  • 

MtjZfj  =  k(A^)A5  +  c(a?)A5  -  -  M^g  + 

m6z6  =  k^^S  +  '  4555  “  M6g  +  v6 


where 


for  l<i<3  A.  =  z  +  i.  sin  9  -  z.  ,  A.  =  z  +  i.0  cos  9  -  z. 

—  -  1  1  11  1  1 

•  •  • 

1  <  i  <  6  A.  =  z  -  i.  sin  9  -  z.  A.  =  z  -  &  9  cos  9  -  z. 

--1  1  1  5  1  i  1 
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Q  =  moment  about  the  center  of  gravity  of  the  main  frame  pro- 

6 

duced  by  horizontal  forces,  Q  =  V  (H.)(8  +  A. ) 

i=l  1  1 

i,  =  distance  from  center  of  gravity  of  main  frame  to  contact 
point  of  i^b  bogie 

*  til 

k(A. )  =  force -deflection  relation  for  1  bogie  suspension 
1  (fig.  C3) 

c(A. )  =  force -velocity  relation  for  l  bogie  suspension 
(fig.  C4) 

til 

V.  =  resultant  vertical  force  of  spring  segments  of  i 
bogie 

til 

H.  =  resultant  horizontal  force  of  spring  segments  of  i 
1  bogie 

8^  =  “  zi  =  rel&bive  displacement  between  adjacent  bogies 

til 

U.  =  spring  constant  for  i  track  spring;  in  this  study, 

1  for  1  <  i  <  5  ui  =  375  lb/in.  ,  =  600  lb/in. 


30.  Observation  of  photographs  of  the  tank  crossing  the  highest 
obstacle  (18  in.)  revealed  that  the  greatest  pitch  angle  expected  would 
be  on  the  order  of  9  deg  or  less.  It  is  seen  that  if 

9=9  deg 

then 

cos  0  =  cos  9  deg  =  0.988  «  1 

and 

9  deg  =  ^  =  0.157  radian 

Since 

sin  9  deg  =  0.156  «  0.157 

the  small  angle  assumption,  i.e.  cos  0  =  1,  sin  0  =  0,  is  valid.  For 
this  reason,  and  to  simplify  the  calculations,  the  small  angle  concept 
was  used  in  the  equations  above. 

31.  Once  the  motions  at  the  center  of  gravity  of  the  main  frame 
have  been  determined,  the  motions  in  the  vicinity  of  the  driver  can  be 
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Fig.  Ch.  m60A1  suspension  damping  force  versus  velocity 

determined  by  combining  the  translatory  and  rotational  motions  to  yield 
the  equation: 


The  value  of  35  in.  represents  the  distance  from  the  center  of  gravity 
to  a  point  2  ft  behind  the  driver,  a  point  at  which  an  accelerometer 
can  be  conveniently  located. 

Computation  of  track  forces 

32.  The  track  compliance  is  represented  chiefly  by  interconnect¬ 
ing  linear  springs  between  the  bogies  and  three  "feelers"  that  are  con¬ 
nected  to  the  front .bogie  by  a  stiff  spring.  The  spring  constants  used 
in  this  study  were  determined  by  observing  photographs  of  the  tank  in 
different  positions  on  an  obstacle.  From  these  photographs,  estimates 
were  made  of  the  influence  on  the  displacement  of  adjacent  bogies  of 
displacing  a  particular  bogie.  Knowing  the  approximate  mass  of  eacn 
bogie  assembly,  an  appropriate  spring  constart  could  be  determined. 

Close  observation  further  revealed  that  upon  approaching  an  obstacle 
larger  than  about  6  in.  high  the  initial  track-obstacle  contact  tended 
to  lift  the  front  bogie  and  guide  it  over  the  obstacle.  This  lift¬ 
ing  has  a  significant  effect  on  the  longitudinal  motion.  To  simulate 
this  effect,  three  feelers  were  positioned  in  front  of  the  first  bogie, 
each  at  a  different  threshold  height,  to  conform  with  the  geometry  of 
the  leading  portion  of  the  track.  The  influence  of  the  feelers  in 
lifting  the  front  bogie  depends  on  the  height  and  shape  of  the  encoun 
tering  obstacle.  Since  no  information  was  available  to  enable  the  de¬ 
termination  of  an  effective  spring  constant,  an  arbitrary  value  of 
600  lb /in.  was  chosen.  Although  this  value  was  estimated,  it  shows 
that  with  a  proper  determination  of  a  spring  constant  these  longitu¬ 
dinal  accelerations  can  be  adequately  simulated. 

Bogie  spring  segments 

33.  The  segmented  wheel  concept  was  used  in  the  model  to  (a> 
enable  predictions  of  longitudinal  accelerations,  (b)  include  important 
geometry  effects  of  the  bogies,  and  (c)  incorporate  a  means  for  de¬ 
scribing  the  composite  compliance  of  the  real  bogie -track -obstacle  sys¬ 
tem.  This  composite  compliance  of  the  real  system  includes  such  phe¬ 
nomena  as  the  effects  of  the  small  terrain  and  obstacle  deformations, 
track  deformations,  and  envelopment  characteristics  among  others  that 
otherwise  in  the  model  would  be  represented  as  infinitely  rigid. 
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34.  Each  bogie  was  divided  into  twelve  10-deg  segments,  six  on 
each  side  of  the  vertical  position  as  shown  in  fig.  C5.  To  account  for 


DIRECTION  OF  TRAVEL 


Fig.  C5.  Schematic  showing  bogie  spring  segment  configuration 


the  track  thickness,  an  effective  bogie  radius  of  17.2  in.  was  used. 
At  these  conditions,  an  average  horizontal  spacing  of  2.357  in.  coin¬ 
cided  closely  with  each  spring  position. 

Determination  of  threshold  heights 

35.  To  compute  the  deflections  of  each  spring  segment,  the  seg 
ment  threshold  heights,  T^  ,  were  first  computed.  These  heights  are 
simply  the  heights  to  each  spring  of  the  undeflected  wheel  (see  fig. 
C5).  1)16  segment  deflection,  ,  is  then  computed  by  the  equation: 

•  zi  ’  h  '  Ti  -  i°| 

0  ’  Yi  -  Ti  -  zi  <  °f 
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where 


Yi  '  vertical  obstacle  height  beneath  the  ith  segment 
\  =  vertical  axle  displacement  of  ith  bogie 

36.  The  segment  deflections  are  permitted  to  have  positi/e  val¬ 
ues  only;  negative  values  are  replaced  by  zero.  The  reference  from 
which  vertical  displacements  are  measured  is  the  point  that  locates  the 
b°g- 5  0x16  when  the  b°Sie  is  imagined  to  be  rigid  and  in  static  equi¬ 
librium.  Static  deviations  from  this  reference  correspond  to  static 
wheel  deflections,  and  superposed  on  these  static  deflections  are  the 
dynamic  obstacle-induced  deflections. 

Computation  of  verti- 
cal  and  horizontal  forces 

37.  The  resultant  vertical  and  horizontal  forces  on  the  first 
bogie  axle  due  to  the  spring  segment  deflections  are  given  by  equations: 

V1  =  jl  fkv  cos  <V5i 
Hi  ’  j^'^h  s:-n  Vh 

where 


kv  and  \  "  segment  spring  constants  for  the  vertical  and  hori¬ 
zontal  modes,  respectively 

=  angle  of  the  ith  segment  from  the  vertical 


=  vertical  deflection  of  the  i 


.th 


segment 


38.  Generally,  ky  and  1^  would  have  the  same  value.  However 
the  increased  stiffness  noted  in  the  horizontal  mode  warranted  a  higher 
spring  constant  in  the  horizontal  mode.  These  values,  k  and  k, 
were  determined  by  examining  oscillograph  records  for  several  tests' 
over  a  10-in. -high  obstacle.  These  k  values  were  adjusted  until  the 
model  outputs  adequately  portrayed  the  gross  features  of  the 
acceleration-time  histories  from  the  oscillographs.  The  choice  of  the 
10-in.  obstacle  was  quite  arbitrary;  it  was  chosen  solely  because  it 
was  close  to  the  median  obstacle  height  that  the  M60A1  tank  would  be 
expected  to  traverse.  Ideally,  such  k  values  would  be  determined  by 
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appropriately  instrumenting  the  bogies  and  performing  a  series  of  sys¬ 
tematic  tests  that  would,  lead  to  the  determination  of  a  proper  set  of 
values . 

39.  Defining  y i  -  ky  cos  0^  and  sin  $  ,  a  y -array 

and  a  a -array  are  established  in  the  same  manner  as  the  array  of 
threshold  heights,  thus  simplifying  the  computations  of  the  resultant 
vertical  and  horizontal  forces. 


Computation  of  moment  pro- 
duced  by  horizontal  forces 

40.  An  important  contribution  to  the  moments  about  the  center  of 
gravity  of  the  main  frame  is  the  horizontal  forces  acting  on  the  bogies. 
A  schematic  diagram,  showing  only  the  front  and  rear  suspensions,  is 
given  in  fig.  C 6  to  illustrate  the  manner  in  which  the  horizontal 
forces  contribute  to  the  moment.  Using  the  small  angle  assumption  and 
the  established  sign  convention  that  the  suspension  deflection,  A  ,  is 
negative  in  the  equilibrium  position,  the  following  equation  is  used  to 
compute  the  moment,  Q  ,  due  to  horizontal  forces. 


q-  -  y  (k.)  t(s 

A  1 


V 


(V»J 


l 


(H.) 


C(s 


U.)e! 


where 

the  quantities  in  brackets  represent  the  moment  arm  frrm  the 

til 

center  of  gravity  to  the  i  bogie 

S  =  vertical  distance  from  the  center  of  gravity  to  the  bogie 
axles  in  the  undeflected  state 

=  resultant  horizontal  force  of  the  l  bogie 

H.  =  longitudinal  distance  from  the  center  of  gravity  to  the 

.  th  ,  . 

l  bogie 

Computer  program 


4l.  A  complete  listing  of  the  digital  computer  program  and  a 
dictionary  of  the  computer  variables  used  are  given  at  the  end  of  this 
appendix.  This  is  a  FORTRAN  IV  program  written  for  a  GE-430  time¬ 
sharing  system.  Mixed  mode  operations,  which  are  acceptable  in  this 
system,  were  used  on  occasion  where  it  proved  advantageous  in  reducing 
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MOMENT  AT  CENTER  OF  GRAVITY  OF  FORCE,  Hj  =  Hj(m  +  n) 

WHERE 

S  +  A  .  ,  . 

n=-E5??  =  s  +  A 

m  =  fj  SIN  9  -  ?,0 
Q,  =  H,[(S  +  A,)  +  ?j0| 

Fig.  C 6.  Schematic  showing  moment  of  horizontal  forces 

the  logic  or  the  number  of  required  statements.  The  file  entitled 
"Fimake"  serves  as  a  convenient  method  to  build  and  input  the  obstacles. 


!£.  The  development  of  the  mathematical  model  of  the  M37  truck 
to  simulate  the  dynamic  response  of  the  truck  while  crossing  rigid  ob¬ 
stacles  followed  the  same  reasoning  as  that  given  in  paragraph  2 6  for 
the  M60A1  tank. 
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43*  A  schematic  diagram  of  the  system  that  was  modeled  is  given 
in  fig.  C7.  This  model  consists  of  four  degrees  of  freedom  that 

Z 


include  bounce  and  pitch  of  the  cr-nter  of  gravity  of  the  main  frame  and 
the  vertical  motions  of  both  axle.:.  A  significant  difference  in  the 
two  mathematical  models  is  that  longitudinal  motion  is  not  taken  into 
account  in  the  M37  truck  model.  The  motions  in  the  vicinity  of  the 
driver*  are  computed  from  the  predicted  motions  at  the  center  of  grav¬ 
ity.  The  frame  of  the  truck  was  considered  rigid,  and  only  the  pneu¬ 
matic  tires  and  suspensions  were  considered  to  contribute  to  the  sprung 
motion  of  the  frame.  The  model  includes  all  pertinent  nonlinearitie^ 
in  the  suspension.  Vehicle  characteristics  used  are  given  in  table  C2. 


*  Accelerometers  can  be  conveniently  located  underneath  the  center  of 
the  driver's  seat. 
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Equations  of  motion 

44.  The  differential  equations  describing  the  motions  of  this 
system  were  developed  in  a  manner  analogous  to  that  described  in  para¬ 
graph  29.  Again,  using  Newton's  second  law  of  motion  and  summing  the 
forces  on  the  body  and  on  each  axle  led  to  the  series  of  equations 
listed  below. 

Forces  on  body  (sprung  mass): 

Mz  =  W^l  -  z  ~  a  sin  0)  +  C1(AI)(z1  -  z  -  ae  cos  e) 

+  ~  z  +  b  sin  0)  +  C2(«^)(z2  -  z  +  b0  cos  0)  -  Mg 

10  -  It^Aj^aCz^  -  z  -  a  sin  0)  +  C^(A^)a(z1  -  z  -  a0  cos  0) 

"  z  +  b  sin  0)  “  C2(^)b(z2  -  z  +  b0  cos  0) 

Forces  on  front  axle  (unsprung  mass): 

Vl  =  'WW  -  z  -  a  sin  0)  -  C1(AL)(z1  -  z  -  a0  cos  0) 

10 

+  £  ’ll**!!  '  "l’  -  "lE 
Forces  on  rear  axle  (unsprung  mass): 

M2Z2  ”  ”^2^^Z2  "  z  +  b  sin  &)  "  c2(/^)(z2  -  z  +  b0  cos  0) 
iO 

*  jT  ri2(pia  -  z2;  - 

where 

~  ?’i  “  7  “  a  s^n  ®  ^  =  z-^  -  7,  -  a0  cos  0 

•  •  •  • 

^3  ~  z2  ”  z  b  s^n  ^  -  Zj  -  z  +  b0  cos  0 

For  this  study,  the  suspension  spring  coefficients  were  represented  by 
third-order  polynomials,  as  shown  in  fig.  c8.  These  polynomials  were 
obtained  by  curve -fitting  the  actual  force -deflection  relations,  and 
are  reasonable  approximations.  The  suspension  damping  is  as  follows: 
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C^(A^)  =  11.8  lb-sec/in. 

=  22.8  lb-sec/in. 
Cg(^)  =  12.0  lb -sec/in. 


(compression) 

(extension) 

(compression) 


=  49*0  lb-sec/in.  (extension) 

Determination  of 
tire -terrain  compliance 

45.  Each  wheel,  which  represented  a  9.00x16,  8-FR  tire  at  45- 
psi*  inflation  pressure,  was  divided  into  ten  segments,  five  on  each 
side  of  the  tire's  center  line,  as  shown  in  fig.  C9.  The  measured 
load-deflection  relation  (fig.  CIO)  for  the  9.00xl6  tire  at  45-psi  in¬ 
flation  pressure  was  such  that  a  center-line  deflection  of  1.5  in.  re¬ 
quired  a  load  of  2860  lb.  At  this  deflection,  four  spring  segments  are 
influenced,  two  on  each  side  of  the  center  line  (fig.  C9).  The  spring 
constant  K  can  be  determined  from  the  statics  equation: 


10 

F  =  Y,  K  cos  0.A. 
i=l  1  1 

where 


\  =  Yi  -  THRESH(i)  -  Z  ,  Y  -  THRESH(i)  -  Z  >  0 

=0  5  Yi  "  THRESH(i)  -  Z  <  0 

Y.^  =  vertical  height  ol  terrain  profile  beneath  ith 
segment 

Z  =  vertical  displacement  of  axle 

THRESH(i)  =  height  from  the  zero  reference  to  the  ith  spring  of 
the  undeflected  wheel  (see  fig.  C9) 


For  this  case  and  due  to  the  symmetry  of  the  segments  about  the  center 
line  the  equation  reduces  to 


*  The  difference  in  the  load -deflection  relation  between  40  psi  (at 
which  subsequent  tests  were  conducted)  and  45  psi  was  deemed 
negligible. 
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2860  =  2K  J  1.4  cos  +  0.7  cos  11.7°  + 


i=l 

where  the  effective  radial  deflections  are: 

=  1.4  in. 
A^  =  A^  =  0.7  in. 


Solving  for  K  yields 


K  =  675  lb/in. 


Defining  GAMMA  =  K  cos  =  675  cos  0^  yields  the  following  relations 
for  the  segments  of  the  front  and  back  wheels: 


GAMMA(l)  =  GAMMA (10) 
GAMMA(2)  =  GAMMA(9) 
GAMMA(3)  =  GAMMA(8) 
GAMMA(4)  =  GAMMA(7) 
GAMMA(5)  =  GAMMA(b) 


-  GAMMA(ll)  =  GAMMA(20)  =  411.75 
=  GAMMA(12)  =  GAMMA(19)  =  506.25 
=  GAMMA(13)  =  GAMMA(18)  =  587.25 
=  GAMKA(14)  =  GAMMA(17)  =  648.00 
=  GAMMA(15)  =  GAMMA(l6)  =  668.25 


A  similar  relation  is  derived  for  the  threshold  heights  of  each  segment 
THRESH(i). 

46.  A  mean  spacing  of  3*07  in.  was  determined  to  be  adequate  for 
portraying  the  projected  spacing  of  the  springs.  As  a  result,  all  pro¬ 
file  points  were  spaced  3.07  in.  apart,  and  no  interpolation  scheme  was 
employed  to  estimate  elevation  between  adjacent  points. 

47.  No  damping  was  incorporated  in  the  tire  compliance  since  in 
actual  vehicles  this  damping  is  negligible  compared  with  that  of  the 
suspensions.  This  truck  model  was  forced  to  traverse  each  obstacle  at 
a  constant  speed,  and  the  outputs  consisted  of  motions  of  the  main 
frame  and  axles  in  terms  of  displacement-,  velocity-,  and  acceleration¬ 
time  histories . 

Gomputer  program 

48.  A  complete  listing  of  the  digital  computer  program  and  a 
dictionary  of  the  computer  variables  used  are  gi\en  at  the  end  of  this 
appendix.  This  is  a  FORTRAN  IV  program  written  for  a  GE-430  time¬ 
sharing  system. 


I 


C26 


t 


PART  III:  TEST  PROGRAMS 

49*  Although  several  series  of  obstacle -vehicle  tests  were  con¬ 
ducted  in  efforts  to  establish  ■‘■esting  procedures  and  criteria  for  dy¬ 
namic  response,  the  first  test  program  that  yielded  data  acceptable  for 
verification  of  the  WES  vehicle  dynamics  model  was  conducted  at  Aber¬ 
deen  Proving  Ground  (APG)  with  an  M60A1  tank.  These  tests  were  con¬ 
ducted  during  the  period  6-9  May  1968.  Subsequently,  an  obstacle  test 
course  was  constructed  at  WES,  and  tests  with  an  M37  truck  were  con¬ 
ducted  intermittently  thereon  during  1968  and  1969  for  another  study. 

Test  Vehicles 

50.  Pertinent  physical  characteristics  for  the  M60A1  tank  and 
the  M37  truck  are  given  in  tables  C3  and  C4,  respectively.  Riotographs 
of  the  vehicles  are  included  in  fig.  Cll.  The  test  vehicles  were 
equipped  with  electronic  systems*  to  measure  and  record  time,  distance 
traveled,  and  vertical  and  longitudinal  accelerations. 

Test  Areas 

APG 

51.  The  test  area  at  APG  was  a  nearly  level  asphalt  strip  (fig. 
C12).  Six  different  sized  obstacles  were  constructed  of  hardwood  and 
fastened  to  the  strip  with  steel  rods.  These  obstacles  were  12  ft  in 
length  and  ranged  from  6  to  18  in.  ii  height.  A  sketch  of  the  obstacle 
configurations  with  dimensions  indicated  is  shown  in  fig.  G13. 

WES 

52.  A  200-ft-long  test  course  (fig.  C14)  constructed  of  concrete, 
asphalt,  and  steel  was  used  for  the  tests  at  WES.  This  course  was  de¬ 
signed  to  accommodate  either  single  or  multiple  obstacles,  both  rigid 
and  deformable.  The  obstacles  can  be  arranged  to  excite  primarily  the 

*  This  instrumentation  is  described  in  detail  in  reference  20. 
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b.  M37  truck,  cargo,  3/4-ton,  4x4 


A<  V?  i  . 


JidfSi 


Fig.  C12.  Test  course  layout,  APG 


pitch  and  bounce  motions  of  the  vehicle  (as  was  done  for  the  tests  re¬ 
ported  herein)  or  to  excite  all  components  of  motion,  i.e.  pitch,  roll, 
bounce,  and  yaw,  in  the  direction  of  and  perpendicular  to  the  direction 
of  travel.  Only  single  obstacles  were  used  in  this  study.  The  ob¬ 
stacles  were  half-round  in  configuration  and  were  welded  directly  to 
steel  rods  embedded  in  the  test  course.  Obstacle  heights  were  6,  8, 

10,  and  12  in. 


Tests  Conducted 


Procedures 

53*  M60A1  tank.  The  vehicle  was  positioned  at  right  angles  to 

and  at  a  sufficient  distance  from  the  obstacle  to  enable  the  desired 
speed  to  be  reached  at  least  5  ft  before  striking  the  obstacle.  The 
test  engineer  instructed  the  driver  to  accelerate  to  a  preselected 
speed  and  try  to  hold  this  speed  constant  while  crossing  the  obstacle. 
Some  tests  were  run  at  as  nearly  the  same  speed  as  possible  over  the 
same  obstacle  height  to  determine  repeatability  of  test  results. 

54.  M37  truck.  The  starting  position  for  each  test  was  suffi¬ 
ciently  far  from  the  obstacle  to  permit  the  driver  adequate  distance  to 
reach  and  maintain  a  constant  vehicle  speed  before  striking  the  ob¬ 
stacle.  Several  methods  of  maintaining  constant  speed  while  crossing 
the  obstacle  were  attempted,  and  the  technique  that  appeared  most  ef¬ 
fective  employed  an  engine  tachometer  calibrated  in  miles  per  hour. 

Some  tests  were  rerun  at  the  approximate  speed  of  the  first  test  to  de¬ 
termine  repeatability  of  test  results. 

Number  of  tests 

55.  The  number  of  tests  conducted  over  selected  obstacle  heights 
is  given  below: 


Obstacle  Height 
_ in. _ 

6 

8 

10 


Number  of  Tests  Conducted 
M60A1  Tank  M37  Truck' 

14  10 

3  11 

4  12 

(Continued) 
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Obstacle  Height 

Number  of  Tests  Conducted 

in. 

M60A1  Tank 

M37  Truck 

12 

11 

2 

16 

9 

* 

18 

_2 

Total 

43 

35 

*  The  vertical  acceleration  was  so  great  on  the 
12-in. -high  obstacle  that  it  was  felt  tests  on 
higher  obstacles  might  damage  the  vehicle. 


Data  Collected 

5 6.  For  each  test,  instrumentation  installed  on  the  test  ve¬ 
hicles  recorded  continuous  measurements  of  time,  distance  traveled  by 
the  vehicle,  and  vertical  and  longitudinal  accelerations  at  selected 
positions  within  the  vehicle.  Event  marks  on  the  oscillogram  indicated 
the  beginning  and  end  of  the  test.  Impact  speed,  peak  vertical  accel¬ 
eration,  and  peak  longitudinal  acceleration  were  determined  from  the 
oscillogram.  A  summary  of  these  data  is  given  in  tables  C5  and  C 6  for 
the  M60A1  tank  and  the  M3 7  truck,  respectively.  Other  data  collected 
included  notes,  observations,  and  photographs. 

Prediction  Data  Obtained 


57.  Predictions  of  peak  vertical  acceleration  and  peak  longitu¬ 
dinal  acceleration  were  made  for  20  obstacle  height -impact  speed  com¬ 
binations  for  the  M60A1  tank,  and  predictions  of  peak  vertical  accel¬ 
eration  were  made  for  14  obstacle  height -impact  speed  combinations  for 
the  M37  truck,  using  the  previously  described  models.  The  speeds  se¬ 
lected  for  the  M60A1  tank  predictions  in  table  C7  were  selected  on  the 
basis  of  preliminary  reduction  of  field  test  data  and  are  generally 
close,  ill  enough  not  identical,  to  the  measured  speeds  for  the  M60A1 
tests  shown  in  table  C5.  The  speeds  for  the  M37  predictions  in  table 
C8  were  arbitrarily  selected  to  yield  dynamic  response  predictions 
above  and  below  the  2,5-g  peak  vertical  acceleration  tolerance  limit. 
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PART  IV:  ANALYSIS  OF  DATA 


Method  of  Analysis  and  Evaluation  of  Predictions 

58.  In  this  study,  vehicle  performance  was  described  in  terms 
of  peak  vertical  acceleration-speed  relations,  peak  longitudinal 
acceleration -speed  relations,  and  speed-obstacle  height  relations  for 
the  2.5-  and  2.0-g  levels  of  peak  vertical  acceleration  and  peak  longi¬ 
tudinal  acceleration,  respectively.  Peak  vertical  acceleration-speed 
relations  and  peak  longitudinal  acceleration -speed  relations  for  sev¬ 
eral  obstacle  heights  were  developed  from  measured  and  predicted  data. 
These  relations  express  the  peak  vertical  acceleration  and  peak  longi¬ 
tudinal  acceleration  to  be  expected- when  a  vehicle  traverses  an  ob¬ 
stacle  of  a  specific  size  at  a  given  speed.  From  these  relations,  the 
speed -ob stacle  height  relations  that  express  the  speed  at  which  the  ve¬ 
hicle  can  traverse  a  given  obstacle  without  exceeding  the  2.5-g  peak 
vertical  acceleration  level  and  the  2.0-g  peak  longitudinal  accelera¬ 
tion  level  were  developed. 

59  The  prediction  accuracy  of  the  models  was  evaluated  by  (a) 
comparing  the  predicted  peak  vertical  accele rat ion- speed  and  peak  longi¬ 
tudinal  acceleration-speed  relations  with  the  measured  peak  vertical 
acceleration- speed  and  peak  longitudinal  acceleration-speed  relations 
and  (b)  comparing  speeds  at  2.5-g  peak  vertical  acceleration  and  2.0-g 
peak  longitudinal  acceleration  levels  developed  from  the  predicted  re¬ 
lations  with  the  speeds  at  2.5-g  peak  vertical  acceleration  and  2.0-g 
peak  longitudinal  acceleration  developed  from  the  measured  data.  The 
quality  of  the  prediction  accuracy  was  evaluated  on  the  basis  of  the  sim¬ 
ple  statistical  parameters,  root  mean  square  (RMS)*  and  percent  error.** 


*  RMS  is  a  measure  of  the  quality  of  a  relation  in  terms  of  the  RMS 
of  the  deviation 


RMS 


where  I !  -  the  sum  of 5  n  =  number  of  deviations  5  and  d  =  deviation, 
i.e.  predicted  minus  measured. 

**  Percent  error  =  £.re^c^e^  ~  measured  y  iqo. 

measured 
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60.  Other  conditions  peculiar  to  a  particular  part  of  the  anal¬ 
ysis  are  discussed  in  the  appropriate  section. 

Peak  Vertical  Acceleration -Obstacle  Height -Speed  Relations 

61.  As  stated,  vertical  acceleration  performance  was  described 
in  terms  of  peak  vertical  acceleration-speed  relations  for  several  ob¬ 
stacle  heights  and  speed-obstacle  height  relations  for  2.5-g  peak  ver¬ 
tical  acceleration  level.  Peak  vertical  acceletation-speed  relations 
developed  from  measured  and  predicted  data  at  the  driver's  seat  are 
shown  in  plates  Cl  and  C2  for  the  M60A1  tank  and  M37  truck,  respec¬ 
tively.  The  curves  drawn  represent  the  lines  of  best  visual  fit. 

Where  the  location  of  the  line  was  doubtful,  judgment  was  aided  by  the 
location  and  curvature  of  lines  on  better  defined  plots.  These  curves 
are  summarized  for  the  M60A1  tank  and  the  M37  truck  in  plates  C3  and 
C4,  respectively.  Speed-obstacle  height  relations  for  both  vehicles 
are  given  in  plate  C5.  The  speed-obstacle  height  relations  were  estab¬ 
lished  from  values  of  speed  and  corresponding  values  of  obstacle  height 
in  plates  C3  and  C4  at  which  2.5-g  peak  vertical  acceleration  occurred. 
M60A1  tank 

62 .  It  can  be  seen  in  plates  Cl  and  C3  that  both  measured  and 
predicted  peak  vertical  accelerations  for  the  M60A1  tank  at  the  driv¬ 
er's  seat  increased  with  an  increase  in  speed.  There  appears  to  be  a 
tendency  for  the  curves  to  crest  at  higher  speeds,  suggesting  that  af¬ 
ter  a  critical  speed  has  been  reached,  a  further  increase  in  speed 
would  not  result  in  an  increase  in  peak  vertical  acceleration.  In 
plate  Cl  it  can  be  seen  that  except  for  the  6-in. -high  obstacles  the 
predicted  peak  vertical  acceleration  was  higher  than  that  measured  at 
low  speeds,  and  lower  than  that  measured  at  higher  speeds.  The  reverse 
is  indicated  for  the  6-in. -high  obstacles.  The  agreement  of  the  meas¬ 
ured  and  predicted  curves  for  the  8-  and  10-in.  obstacles  is  very  good, 
and  at  the  2.5-g  tolerance  limit  the  curves  for  the  12-,  16-,  and  l8-in. 
obstacles  seem  to  be  reasonably  close. 

63.  In  plate  C5,  the  speed-obstacle  height  curve  for  the  M6OAI 
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tank  shews  that  the  speed  at  which  2.5-g  peak  vertical  acceleration  is 
reached  at  the  driver's  seat  decreases  with  an  increase  in  obstacle 
height.  The  effect  of  obstacle  height  on  2.5-g  peak  vertical  accelera¬ 
tion  begins  to  diminish  rapidly  at  about  the  9-in.  obstacle  height.  An 
examination  of  predicted  data  points  shows  that  they  are  in  good  agree- 
ment  with  curves  developed  from  the  measured  data. 

M37  truck 

64.  In  plates  C2  and  C4,  both  measured  and  predicted  peak  verti¬ 
cal  accelerations  for  the  M37  truck  at  the  driver’s  seat  increased  with 

“  increase  in  sPeed-  At  the  higher  speeds  the  rate  of  change  in  peak 
vertical  acceleration  was  less  than  that  at  lower  speeds.  Except  in  a 
few  cases  at  very  low  speeds  (plate  C2),  there  is  a  large  difference 
between  the  measured  and  predicted  peak  vertical  accelerations  at  the 
same  speed,  with  the  predicted  peak  vertical  acceleration  being  lower 
than  the  measured.  An  explanation  for  this  is  not  readily  available 
without  further  study.  It  may  be  because  of  the  difference  in  the  phys¬ 
ical  characteristics  of  the  vehicle  and  the  model,  e.g.,  spring  rates, 
damping  rates,  etc.,  or  because  of  the  inability  to  accurately  produce 
a  computer  simulation  of  the  motion  history  of  the  vehicle  immediately 
prior  to  and  during  obstacle  traversal. 

65.  The  speed-obstacle  height  curve  for  the  M37  truck  in  plate 
C 5  was  established  from  the  measured  relations  for  the  6-,  8-,  and 
10-in. -high  obstacles.  (Note  in  plate  C2  that  even  at  low  speeds  the 
peak  vertical  acceleration  measured  over  12 -in. -high  obstacles  greatly 
exceeded  2.5  g’s.)  The  relation  is  similar  to  that  discussed  above  for 
the  M60A1  tank  except  that  2.5-g  peak  vertical  acceleration  at  the 
driver's  seat  was  reached  at  lower  speeds  and  obstacle  heights  than 
that  for  the  M60A1  tank.  It  can  be  seen  that  the  predictions  indicate 
a  higher  speed  for  crossing  8-  and  10-in. -high  obstacles  than  that 
shown  by  the  relations  developed  from  measured  data.  In  addition,  the 
predicted  speed  for  the  12-in. -high  obstacle  was  4  mph,  while  the  meas¬ 
ured  data  indicate  (see  again  plate  C2)  that  the  vehicle  cannot  cross  a 
12-in. -high  obstacle  at  any  speed  without  exceeding  2.5-g  peak  vertical 
acceleration.  No  predicted  value  is  shown  for  6-in. -high  obstacles 


C34 


because  the  predictions  of  peak  vertical  acceleration  at  the  driver's 
seat  did  not  reach  2.5  g  at  speeds  up  to  20  mph.  Hence,  it  is  apparent 
that  for  all  obstacle  heights  tested,  the  predicted  speed  at  which  2.5-g 
peak  vertical  acceleration  at  the  driver's  seat  would  occur  was  greater 
than  that  indicated  by  the  relations  developed  from  measured  data. 
Prediction  accuracy 

66.  The  measured  and  predicted  peak  vertical  accelerations,  the 
deviations,  and  the  percent  error  for  the  speeds  at  which  comparisons 
were  made  are  given  in  tables  C7  and  C8  for  the  M60A1  tank  and  the  M37 
truck,  respectively.  The  measured  data  shown  in  column  4  of  tables  C7 
and  C8  were  obtained  from  the  relations  established  for  tfleasured  data 
shown  in  fig.  a,  plate  C3,  for  the  M60A1  tank  and  in  fig.  a,  plate  C4, 
for  the  M37  truck.  The  predicted  data  given  in  column  3  of  tables  C7 
and  C8  are  for  predictions  made  for  the  speeds  shown  in  -olumn  2. 

67.  The  measured  and  predicted  peak  vertical  acceleration  data 
given  in  tables  C7  and  C8  are  plotted  in  plate  C 6  for  both  vehicles, 
and  a  summary  of  the  RMS  and  average  percent  error  for  each  obstacle 
height  is  given  in  the  following  tabulation: 


M60A1  Tank 

N37  Truck 

Obstacle 

Average 

Average 

Height 

Percent 

Percent 

in. 

RMS 

Error 

RMS 

Error 

6 

0.59 

81 

2.06 

65 

8 

0.10 

11 

1.45 

28 

10 

0.22 

41 

1.17 

18 

12 

0.43 

18 

4.00 

68 

16 

0.87 

22 

mi  — 

18 

0.16 

6 

— 

— 

Overall 

0.51 

37 

1.94 

39 

In  the  tabulation  above  it  can  be  seen  that  the  overall  percent  error 
for  the  tracked  vehicle  was  only  slightly  lower  than  for  the  wheeled 
vehicle,  but  in  terms  of  RMS  the  predictions  for  the  tracked  vehicle 
were  much  better  than  for  the  wheeled  vehicle.  For  each  vehicle  at 
several  obstacle  heights  the  percent  error  and  RMS  are  sufficiently 
large  to  indicate  that  some  refinement  is  needed  in  the  mathematical 
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models.  In  fig.  a,  plate  C6,  the  predictions  for  the  M60A1  tank  are 
generally  reasonable  except  for  a  few  points--the  two  highest  accelera¬ 
tions  for  the  6-  and  l6-in.  obstacles,  and  the  highest  acceleration 
for  the  12 -in.  obstacle.  In  fig.  b,  plate  C 6,  the  predictions  for  the 
M37  truck  are  generally  poor  except  for  three  points— the  two  lowest 
accelerations  for  the  10-in.  obstacle  and  the  lowest  acceleration  for 
the  8 -in.  obstacle. 

68.  A  comparison  of  measured  and  predicted  speeds  at  2.5-g  P-ak 
vertical  acceleration  was  made  by  selecting  the  speed  and  obstacle 
heights  at  which  2.5-g  peak  v-rtical  acceleration  occurred  in  the  fig¬ 
ures  given  in  plates  C3  and  C4.  The  data  obtained  in  this  manner, 
along  with  values  of  prediction  accuracy,  are  given  in  the  following 
tabulation. 


Obstacle 

Height 

in. 

M60A1  Tank 

M37 

Truck 

Speed,  i 
2.5-1 

jiph ,  at 
v  PVA* 

Deviation 

( Predicted  Percent 
-  Measured)  Error 

Speed,  mph,  at 
2.5-g  PVA 

Deviation 

( Predicted  Pe  rcent 

-  .assured)  Error 

Predicted 

Measured 

Predicted 

6 

_  _ 

-- 

— 

8.8 

— 

— 

— 

8 

15.6 

15.6 

O.o 

0 

5.1 

7.3 

2.2 

35 

to 

11.0 

u.6 

0.6 

6 

h.b 

5.0 

0.6 

1** 

12 

7.U 

7.2 

-0.2 

3 

— 

u.o 

16 

5-3 

6.1 

0.8 

15 

— 

-- 

18 

It  .9 

5.1 

0.2 

U 

-- 

"" 

" 

Average 

7 

Ave  rage 

25 

*  I'VA  designates  peak  vertical  acceleration. 


The  speed-obstacle  height  relations  developed  from  the  data  in  the  tab¬ 
ulation  above  are  shown  in  plate  C5.  It  can  be  seen  that  the  relations 
are  reasonably  well  defined  and  that  as  obstacle  height  increases,  the 
speed  at  which  the  limiting  acceleration  value  is  reached  decreases  and 
that  the  predicted  data  agree  reasonably  well  with  the  curves  drawn  for 
the  measured  data,  particularly  for  the  M60A1  tank.  For  a  given  ob¬ 
stacle  height,  the  M60A1  tank  can  cross  that  obstacle  at  almost  twice 
the  speed  that  the  M37  truck  can  cross  it  before  2.5-g  peak  vertical  ac¬ 
celeration  is  reached.  The  tabulation  above  and  plate  C7  reveal  that 
the  accuracy  of  predicting  speeds  at  2.5-g  vertical  accleration  is  very 
good  for  the  M60A1  tank  and  reasonably  good  for  the  K37  truck  throughout 

the  speed  range  shown. 
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^gak  Longitudinal  Acceleration -Obstacle  Height -Speed  Relations 


69.  Longitudinal  acceleration  performance  was  described  in  terms 
of  peak  longitudinal  acceleration-speed  relations  for  several  obstacle 
heights,  and  speed -obstacle  height  relations  for  2.0-g  peak  longitu¬ 
dinal  acceleration.  Peak  longitudinal  acceleration-speed  relations  de¬ 
veloped  from  measured  and  predicted  data  at  the  driver's  seat  are  shown 
in  plate  C8  for  the  M60A1  tank.  Since  no  predictions  of  peak  longitu¬ 
dinal  acceleration  were  made  for  the  M37  truck,  the  peak  longitudinal 
acceleration-speed  relations  from  measured  data  only  are  shown  in  plate 
C9.  The  curves  drawn  represent  the  lines  of  best  visual  fit,  and  they 
are  summarized  in  plates  C9  and  CIO  for  the  M37  truck  and  the  M60A1 
tank,  respectively.  Speed-obstacle  height  relations  for  the  M60A1  tank 
are  shown  in  plate  Cll.  The  speed-obstacle  height  relation  for  the 
M60A1  tank  was  established  from  values  of  speed  and  corresponding  val¬ 
ues  of  obstacle  height  in  plate  CIO  at  which  2.0-g  peak  longitudinal 
acceleration  occurred. 

M60AI  tank 

70.  It  can  be  seen  in  plates  C8  and  CIO  that  both  measured  and 
predicted  peak  longitudinal  accelerations  for  the  M60A1  tank  at  the 
driver's  seat  increased  with  an  increase  in  speed.  In  all  cases,  the 
predicted  curves  indicate  higher  values  of  peak  longitudinal  accelera¬ 
tion  than  the  measured  curves.  Although  the  scatter  in  data  is  not 
excessive,  additional  data  are  required  to  better  define  the  relations. 
In  table  C5  it  can  be  seen  that  the  cest  data  include  only  one  measured 
peak  .longitudinal  acceleration  greater  than  2.0  g's.  However,  the 
curves  for  the  larger  obstacles  shown  in  plates  C8  and  CIO  were  ex¬ 
tended  to  2.0  g's. 

71.  The  speed-obstacle  height  relations  (plate  Cll)  for  the 
M60A1  tank  show  that  the  speed  at  which  2.0-g  peak  longitudinal  accel¬ 
eration  is  reached  at  the  driver's  seat  decreases  with  an  increase  in 
obstacle  height.  The  predicted  data  shown  in  the  plot  do  not  agree 
very  well  with  the  measured  data. 
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M37  truck 


72.  The  curves  in  plate  C9,  developed  for  measured  data  only,  in¬ 
dicate  that  the  peak  longitudinal  acceleration  of  the  M37  truck  increased 
with  speed  up  to  about  8  to  10  mph,  and  then  decreased  with  further  in¬ 
crease  in  speed.  An  increase  in  peak  longitudinal  acceleration  with 
increased  obstacle  height  is  also  indicated,  as  would  be  expected. 

While  the  scatter  of  the  data  points  is  not  excessive,  there  is  an 
obvious  need  for  additional  data,  especially  for  the  12 -in.  obstacle 
height . 

Prediction  accuracy 

73.  The  quality  of  the  prediction  accuracy  for  peak  longitudinal 
acceleration-speed  and  obstacle  height-speed  relations  for  the  M60A1 
tank  was  determined  in  the  same  manner  as  that  for  the  vertical  accel¬ 
eration  relations.  The  measured  and  predicted  data,  deviation,  and  per¬ 
cent  errors  for  the  speeds  at  which  comparisons  were  made  are  given  in 
table  C7.  The  measured  data  shown  in  column  8  of  table  C7  were  ob¬ 
tained  from  the  relations  established  from  measured  data  shown  in  plate 
CIO,  The  predicted  data  shown  in  column  7  are  for  predictions  made  for 
the  speeds  shown  in  column  2. 

74.  The  measured  and  predicted  peak  longitudinal  acceleration 
data  given  in  table  C7  for  the  M60A1  tank  are  plotted  in  plate  C12,  and 
a  summary  of  the  RM/>  and  average  percent  error  for  each  obstacle  height 
is  given  below. 


M6OAI  Tank 


Obstacle 

Height 

in. 

RMS 

Average 

Percent 

Error 

6 

0.26 

74 

8 

0.33 

80 

10 

0.22 

43 

12 

0.56 

102 

16 

1.07 

119 

18 

1.05 

76 

Overall 

0.64 

82 
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The  preceding  tabulation  shows  that  the  overall  prediction  accuracy  is 
not  very  good  and  no  consistent  pattern  is  evident  other  than  that  pre¬ 
diction  accuracy  generally  decreased  with  obstacle  height  increase.  In 
plate  C12  the  individual  data  points  reveal  perhaps  more  meaningful 
trends.  Except  for  two  points  for  the  6-in.  obstacle  height,  predicted 
peak  longitudinal  accelerations  are  higher  than  the  measured  accelera¬ 
tions.  Predictions  for  the  12-,  16-,  and  l8-in.-high  obstacles  are  con¬ 
sistently  higher  than  the  measured  data.  The  deviations  for  these  ob¬ 
stacle  heights  are  also  much  greater  than  for  the  other  obstacle 
heights . 

75.  A  comparison  of  M60A1  measured  and  predicted  speeds  at  2.0-g 
peak  longitudinal  acceleration  was  made  by  selecting  the  speed  and  ob¬ 
stacle  height  at  which  2.0-g  peak  longitudinal  acceleration  occurred  in 
plate  CIO.  Since  only  two  points  were  available  for  the  2.0-g  level  of 
acceleration,  a  comparison  was  also  made  for  the  1.0-g  level  of  peak 
longitudinal  acceleration.  The  data  obtained  in  this  manner,  along 
with  values  of  prediction  rc curacy,  are  given  below. 


M6OAI  Tank 


Obstacle 

Height 

in. 

°peed,  mph,  at 
2.0-g  FLA* 

Deviation 
(Predicted 
-  Measured) 

Percent 

Error 

Speed,  mph,  at 
1.0-g  HA 

Deviation 
(Predicted 
-  Measured) 

Percent 

Error 

Measured 

Predicted 

Measured 

Predicted 

6 

-- 

— 

-- 

— 

17.3 

... 

-- 

— 

8 

-- 

-- 

-- 

-- 

U.  5 

12.0 

-2.5 

-17 

10 

-- 

— 

-- 

— 

12.0 

10.6 

-1.1* 

-12 

12 

-- 

— 

— 

— 

8.0 

I*. 2 

-3.8 

-1*8 

16 

10.0 

5.6 

-!*.!* 

-1*1* 

6.6 

-- 

-- 

-- 

18 

7.U 

it.6 

-2.8 

-38 

!*.6 

— 

— 

— 

Average 

-1*1 

Average 

-26 

*  PIA  designates  peak  longitudinal  acceleration. 


The  speed-cV Stacie  height  relations  developed  from  the  data  given  in  the 
tabulation  above  for  the  M60A1  tank  are  shown  in  plate  CU.  The  data 
from  which  the  curve  for  the  speed-obstacle  heights  at  which  2.0-g  peak 
longitudinal  acceleration  is  reached  are  limited,  but  the  curve  for  the 
1.0-g  level  is  well  defined.  Both  curves  show  that  as  the  obstacle 
height  increases,  the  speed  at  which  the  limiting  acceleration  value 
is  reached  decreases.  From  the  data  above  and  plate  C13,  it  can  be 
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seen  that  the  predicted  speed  at  which  a  given  level  of  peak  longitu¬ 
dinal  acceleration  will  occur  for  a  given  obstacle  height  is  lower  than 
the  measured  speed.  The  difference  in  measured  and  predicted  values  is 
greater  for  the  2.0-g  level  than  for  the  1.0 -g  peak  longitudinal  accel¬ 
eration  level.  The  prediction  accuracy  for  the  1.0-g  peak  longitudinal 
acceleration  level  is  acceptable  for  the  8-  and  10-in.  obstacle  heights, 
but  it  is  not  acceptable  for  the  12 -in.  obstacle.  For  the  2.0-g  peak 
longitudinal  acceleration  level,  the  prediction  accuracy  is  not 
acceptable. 


Notes  and  Observations 


7 6.  During  the  test  program,  it  was  observed  that  in  all  of  the 
tests  the  vehicles  appeared  to  strike  the  obstacles  at  or  very  near  an 
angle  of  90  deg;  however,  in  some  tests  there  was  a  change  in  vehicle 
orientation  during  traversal  of  the  obstacle.,  Since  the  peak  vertical 
acceleration  generally  occurs  when  the  vehicle  strikes  the  ground  after 
crossing  the  obstacle,  the  magnitude  of  the  peak  vertical  acceleration 
is  influenced  by  the  attitude  and  orientation  of  the  vehicle.  The 
M60A1,  for  instance,  might  strike  the  ground  with  both  tracks  simulta¬ 
neously,  thus  deriving  the  maximum  benefit  of  the  suspension  system; 
one  track  may  strike  the  ground  first  with  the  vehicle  in  such  position 
that  the  major  portion  of  the  shock  is  imposed  on  the  suspension  system 
of  one  track  only,  resulting  in  some  roll  motion;  or  the  vehicle  may 
strike  the  ground  in  any  position  between  these  two  extremes  with  a 
possibility  of  a  slight  change  in  vehicle  direction.  The  extremes  are 
even  wider  for  the  M37  truck.  A  wheeled  vehicle  may  land  on  all  wheels 
simultaneously,  on  two  wheels,  or  on  a  single  wheel.  Not  all  of  these 
extremes  were  evidenced  in  these  tests;  however,  there  were  occasions 
when  one  tr'  :h  or  wheel  appeared  to  make  contact  with  the  ground  before 
the  other  track  or  wheel. 

77.  In  the  M60A1  tests,  the  peak  longitudinal  acceleration  at 
the  driver's  seat  occurred  when  the  vehicle  contacted  the  obstacle, 
whereas  the  peak  longitudinal  acceleration  of  the  M37  truck  generally 


occurred  when  the  front  wheels  struck  the  obstacle,  although  there  were 
some  occasions  when  the  peak  longitudinal  acceleration  occurred  when 
the  rear  wheels  contacted  the  obstacle.  Since  some  combinations  of  ob¬ 
stacle  height  and  speed  caused  the  M37  truck  to  become  airborne,  i.e. 
lose  contact  with  the  ground,  it  is  easy  to  hypothesize  a  condition  in 
which  the  peak  longitudinal  acceleration  for  the  JQ7  might  occur  when 
the  vehicle  strikes  the  ground  after  crossing  the  obstacle.  The  M60A1 
tank  did  not  become  airborne  during  any  of  the  tests  described  herein. 

78.  As  previously  stated,  the  drivers  attempted  to  maintain  a 
constant  speed  across  the  obstacle.  Since  this  would  be  patently  im¬ 
possible  if  the  driver  waited  until  he  felt  the  vehicle  slowing  before 
he  applied  additional  power,  the  procedure  evolved  was  to  apply  addi¬ 
tional  power  as  close  to  the  moment  of  impact  as  possible.  The  success 
of  this  procedure,  in  terms  of  repeatability  of  the  test  results,  has 
its  limitations.  Obviously,  there  was  nothing  the  driver  could  do  to¬ 
ward  maintaining  a  constant  speed  when  the  M37  truck  was  airborne. 

Rather  than  attempting  to  maintain  a  constant  speed,  it  is  suggested 
that  future  tests  be  conducted  with  the  throttle  in  a  fixed  position 
throughout  the  test. 
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PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


79*  Based  on  the  data  reported  herein  tmd  subject  to  the  limits 
imposed  by  these  data,  the  following  conclusions  are  offered: 

a.  Discrete  rigid  obstacles  affect  vehicle  performance  by 
producing  adverse  vertical  and  longitudinal  motions 
that  may  endanger  the  driver  or  damage  the  cargo. 

b.  The  magnitudes  of  peak  vertical  and  longitudinal  accel¬ 
erations  are  dependent  primarily  on  speed,  obstacle 
height,  and  characteristics  of  the  vehicle  suspension 
system. 

c_.  Vehicle  performance  in  terms  of  peak  vertical  and  peak 
longitudinal  accelerations  when  traversing  discrete, 
rigid  obstacles  can  be  correlated  with  impact  speed. 

d.  The  speed  at  which  the  M60A1  tank  may  contact  a  rigid 
obstacle  without  exceeding  the  2.5-g  peak  vertical  ac¬ 
celeration  tolerance  limit  or  the  2.0-g  peak  longitu¬ 
dinal  acceleration  tolerance  limit  can  be  predicted  by 
the  described  mathematical  techniques  with  reasonable 
accuracy. 

e.  The  speed  at  which  the  M37  may  contact  a  rigid  obstacle 
without  exceeding  the  2.5-g  peak  vertical  acceleration 
tolerance  limit  can  be  predicted  by  the  described  mathe¬ 
matical  techniques  with  reasonable  accuracy  for  8-  to 
10-in. -high  obstacles.  Above  and  below  this  height 
range,  however,  the  prediction  accuracy  left  much  to  be 
desired. 

f.  Refinement  is  needed  in  the  dynamic  response  prediction 
models  with  more  improvement  required  in  predicting  dy¬ 
namic  response  of  wheeled  vehicles. 


Recommendations 


80.  It  is  recommended  that: 

a.  Additional  systematic  controlled  testing  be  done  with 
tracked  and  wheeled  vehicles  to  refine  and  extend  the 
relations  and  revise  the  vehicle  dynamic  prediction 
models  presented  herein. 

b_.  Refinements  be  made  in  the  dynamic  response  prediction 
models  to  bring  the  predictions,  particularly  those  for 


higher  speed  and  acceleration  levels,  more  clearly  in 
line  with  the  measured  dynamic  response. 

Investigation  be  continued  to  explore,  in  addition  to 
peak  acceleration,  other  performance  parameters  such  as 
root  mean  square  acceleration  and  absorbed  power  as  tol¬ 
erance  descriptors. 
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Table  Cl 


Vehicle  Characteristics  Used  in  Predicting  m60A1  Dynamic  Response 


Distance  from  road  wheel  to  body  center  of  gravity: 

Road  wheel  . 

Road  wheel  . 

Road  wheel  . 

Road  wheel  . 

Road  wheel  i . 

Road  wheel  i g . . 

Body  pitch  inertia  . 

Sprung  weight  . 

Vertical  distance  S  from  each  bogie  center 

to  body  center  of  gravity  in  no-load  condition  . 

Effective  unsprung  weight  for  each  bogie, 

lower  track  assembly  . 


.  7 6  in. 

.  44  in. 

.  12  in. 

.  24  in. 

.  56  in. 

.  88  in. 

581,700  in. -lb/sec2* 
.  .  .  .  54,500  lb* 

.  46  in. 

.  1420  lb 


*  The  values  for  the  body  pitch  inertia  and  sprung  weight  represent  but 
one-half  of  the  true  value  in  order  to  fit  the  two-dimensional  model 
which  represents  one-half  of  the  vehicle  that  is  assumed  to  be  split 
down  its  longitudinal  axis. 


Table  C2 


Vehicle  Characteristics  Used  in  the  Dynamic 
Response  Predictions  for  the  M37  Truck 

Vertical  distance  from  wheel  hub  center  to  body  center  of  gravity: 


Road  wheel  .  20.95  in. 

Road  wheel  &2 .  21.50  in. 

Road  wheel  ^ .  21.00  in. 

Road  wheel  . - .  21.95  in. 

Undeflected  wheel  radius . .  .  17.0  in 

Total  weight  (loaded)  .  7550  lb 

Sprung  weight .  6212  lb 

Unsprung  weight .  133q  lb 

Front  axle .  722  lb 

Rear  axle .  6l6  lb 

Sprung*  pitch  moment  of  inertia  about 


center  cf  gravity  .  50,892  ib-sec2/in. 

Sprung  roll  moment  of  inertia  about 

center  of  gravity  .  11,425  ib-sec2/in. 

Front  axle  roll  moment  of  inertia .  1256  lb-sec2/in. 

Rear  axle  roll  moment  of  inertia .  955  lb-sec2/in. 

Wheel  travel  from  static  to  bump  stop: 

Front .  4  in. 

Rear .  8.5  in. 

Tire  (9. 00-16,  8 -PR,  45  psi)  damping  rate .  9.56  lb-sec/in. 


Note:  One-half  of  these  values  was  used  in  the  two-dimensional  model 
under  the  assumption  that  the  vehicle  is  symmetrical. 

*  Sprung  mass,  M  =  =  16.1  lb-sec2/in. 

Front  axle  unsprung  mass,  Mx  =  |||  =  I.87  lb-sec2/in. 

Rear  axle  unsprung  mass,  M2  =  fgf  =  1-60  lb-sec2/in. 


Table  C 3 

Summary  of  Vehicle  Characteristics  and 
Performance  Data  for  M6OAI  Tank 


Cross-country  gross  weight,  fully  equipped  plus  payload  and 
personnel 


Track  weight,  lb 

Left . 

Right . 

Total  gross  weight,  lb 
Payload,  lb . 

Dimensions 


•  •  .  .  54,500 

•  •  •  •  54,500 
.  .  .  .  109,000 
6,000  plus  crew' 


Overall  length,  in.  (hull  only) . 

Height  of  leading  edge,  in . 

Width,  in .  [ 

Horizontal  distance  from  leading  edge  (omitting  gun)  to 
center  line  of  front  sprocket  of  vehicle,  in.  ... 

Vehicle  approach  angle,  deg  . 

Vehicle  departure  angle,  deg  . 

Ground  clearance  of  hull  between  tracks,  in . 

Force  leading  edge  can  withstand,  lb . . 

Winch  capacity,  lb  . 

Water  performance  characteristics 


273.5 

45.0 

143.0 

21.77 

90 

90 

18 

limited 

None 


Fording  depth 

Normal  fording  (no  kit),  in.  . 

Fording  with  kit,  in . 

Fording  with  snorkel  kit,  in. 

Water  speed,  mph  . 

Volume  of  water  displaced,  cu  ft 

Engine 


48 

69 

162 

2 

910.4 


Make . 

Model . 

Fuel  type . 

Brake  horsepower  .  . 
Maximum  torque,  lb-ft 
RIM  at  maximum  torque 

Transmission 


.  Continental 
AVDS-1790-2A 
.  .  Diesel 

.  .  642 

•  •  1720 

.  .  1800 


Make  .  .  . 
Model.  .  . 

Gear  ratios 
1st  gear 


.  Allison 
CD-850- 6A 

.  3.497:1 


(Continued ) 


Table  C3  (Concluded) 


2nd  gear  . 

Torque  convertor  ratio  (maximum) 
Final  drive 

Type . 

Ratio  . 

Pitch  diameter  of  sprocket,  in. 

Track  data 


.  1.256:1 

•  •  •  3.43:1  (computed) 

Single  stage  speed  reduced 

.  5.08:1 

.  24.504 


Length  of  track  in  contact  with  ground,  in.  ... 

Width  of  track,  in . 

Grouser  height,  in . [  * 

Total  number  of  bogies  in  contact  with  ground.  .  . 

Area  of  one  track  shoe,  sq  in . 

No.  of  tracks . . 

•••••••••• 

Steering  data 

Turning  radius  at  a  forward  speed  determined  from 
point  at  center  of  outside  track,  ft 

Pivot  . 

Normal  . 

Horizontal  distance  from  leading  edge  of  vehicle  to 
midpoint  of  track  in  contact  with  the  ground,  in. 
Width  of  track,  in . 

Center  of  gravity  location 

Horizontal  distance  from  front  sprocket,  in.  .  . 
Vertical  distance  above  bottom  of  track,  in.  .  . 

Maximum  performance  data 


166.7 

28.0 

1.5 

12 

194.3 

2 


9.6' 

35.0 

135.41 
28. C 


73.45 

54.5 


Maximum  speed,  mph 

Forward . 

Reverse . 

Gradeability,  %  . 

Bridging,  in . * 

Side-slope  holding  capacity,  %  ....... 

Maximum  vertical  height  vehicle  will  climb,  in 


.  .  30.0 

Not  applicable 

.  .  60 

.  .  102 

•  .  36 

.  •  36 


Table  C4 

Summary  of  Vehicle  Characteristics  and 
Performance  Data  for  the  M37  Truck 


Cross-country  gross  weight,  fully  equipped  plus  payload  and 
personnel 

Axle  loads  (front  to  rear),  lb 

No.  1 . 

No.  2 . 

Total  gross  weight,  lb . * 

Payload,  lb  (cross  country)  . 

Dimensions,  in. 

Overall  length  (including  winch  if  available)  .... 

Wheelbase  . 

Height  of  leading  edge . 

Width . *  [  [  ‘ 

Distance  between  wheel  center  lines  . 

Distance  between  axles . 

Vehicle  approach  angle,  deg  . 

Vehicle  departure  angle,  deg  . 

Undercarriage  clearance,  in. 

Axle . 

Interior . . 

•••••••••••• 

Force  leading  edge  can  withstand,  lb  . 

Winch  capacity,  lb  . 

Water  performance  characteristics 

Fording  depth,  in. 

Normal . 

With  fording  kit  . 

Water  speed,  mph . 

Engine 


3,600 

3,600 

7,200 

1,500 


189.4 

112 

26 

72.8 

62 

112 

38 

32 


11 

16 


10,000 

7,50° 


30 

79 

2 


Make . 

Model  . 

Fuel  type  . 

Brake  horsepower . 

Maximum  torque,  lb-ft  .  .  .  .  . 
Engine  RPM  at  maximum  torque 
Engine  RPM  at  brake  horsepower. 


Dodge 

T-245 

Gasoline 

77 

188 

1,200 

3,200 


(Continued) 


Table  C4  (Concluded) 


Transmission 

Type  or  model  . 

Ratios 

1st . 

2d . 

3d . 

4th . . . 

Transfer  case 

Model  . 

Ratios 

High . 

Low . 

Axles 

Model  . 

Ratio . ' . 

Tire  data 

Type . 

Size . 

Ply . 

Tread  design . 

Unloaded  diameter  (including  tread),  in. 

Unloaded  width,  in . 

Tread  depth,  in . 

No.  of  tires . 

No.  of  wheels  . 

Cross-country  inflation  pressure,  psi  .  .  . 
Highway  inflation  pressure,  psi  . 

Steering  data 

Turning  radius  (curb  to  curb),  ft  . 

Time  required  to  steer  from  straight 

ahead  to  full  lock  turn,  sec  . 

Distance  from  front  wheel  steering  (hub) 

pivot  to  front  of  vehicle,  in . 

Distance  from  front  wheel  steering  (hub) 
pivot  to  outside  of  vehicle,  in . 

Center  of  gravity  location,  in. 

Horizontal  distance  from  front  axle 
Vertical  distance  above  ground  at  full  load 
static  position 


T-245-395!? 

6.40:1 

3.09:1 

1.69:1 

1.00:1 


Timken 

1.00:1 

1.96:1 


T-245 

5.83:1 


Mud  and  snow 
9.00x16 
8 

NDMS 

05 

10.2 

1 

4 

4 

14 

4o 


25 

3 

35-62 

10.0 


63 

34 


Table  C5 


Summary  of  Measured  Data  and  Test  Results  for 


M60A1  Tank  Tests 


Field 

Obstacle 

Impact 

Identification 

Height 

Speed 

No. 

in. 

mph 

1-1 

6 

3.75 

1-2 

6 

4.87 

1-3 

6 

3.52 

2-1 

6 

6.20 

2-2 

6 

5.88 

2-3 

6 

5.81 

3-1 

6 

10.51 

3-2 

6 

10.27 

3-3 

6 

10.05 

4-1 

6 

16.50 

4-2 

6 

17.17 

4-3 

6 

17.47 

5-2 

6 

18.15 

5-3 

6 

18.15 

15-1 

8 

6.15 

15-2 

8 

11.25 

15-3 

8 

17.62 

16-1 

10 

6.07 

16-2 

10 

IO.65 

16-3 

10 

17.17 

16-4 

10 

17.27 

6-1 

12 

4.42 

6-2 

12 

5.02 

6-3 

12 

4.87 

7-1 

12 

6.45 

7-2 

12 

5.85 

7-3 

12 

6.07 

9-1 

12 

8.44 

9-2 

12 

8.60 

9-3 

12 

8.60 

8-1 

12 

10.17 

8-2 

12 

IO.52 

10-1 

16 

5.40 

10-2 

16 

4.50 

10-3 

16 

4.42 

11-1 

16 

6. 37 

11-2 

16 

6.91 

11-3 

16 

6.52 

12-1 

16 

8.?, 

12-2 

16 

7.65 

12-3 

16 

8.77 

13-2 

18 

5.40 

13-3 

18 

4.92 

Peak  Vertical 
Acceleration 


Peak  Longitudinal 
Acceleration 


Table  C 6 


44 

F68-0024 

6 

4.0 

1.48 

0.74 

**5 

F68-0025 

6 

3.5 

1.43 

O.76 

46 

F68-0026 

6 

6.4 

1.48 

0.82 

47 

F68-0027 

6 

8.9 

2.57 

NA 

48 

F68-0028 

6 

11.8 

3.40 

0.80 

49 

F68-0029 

6 

11.6 

3.22 

0.76 

50 

F68-0030 

6 

15.1 

3.66 

0.79 

51 

F68-0031 

6 

14.5 

3.40 

0.82 

52 

F68-0032 

6 

15.3 

3.85 

0.70 

53 

F68-0033 

6 

14.8 

3.64 

0.81 

54 

F68-0103 

8 

2.8 

1.64 

0.85 

55 

F68-0104 

8 

3.0 

1.63 

0.96 

56 

F68-0105 

8 

8.1 

4.38 

1.57 

57 

F68-0106 

8 

5-4 

2.43 

1.50 

58 

F68-0107 

8 

10.8 

5.78 

1.48 

59 

F68-0109 

8 

13.3 

6.59 

1.45 

60 

F68-0110 

8 

13.3 

6.20 

1.50 

61 

F68-0111 

8 

19-1 

7-99 

1.19 

62 

F68-0112 

8 

9-3 

3.44 

1.48 

63 

F68-0113 

8 

19.O 

8.73 

1.25 

64 

F68-0114 

8 

19.6 

8.72 

1.4o 

65 

F68-0034 

10 

2.9 

1.68 

I.07 

66 

F68-0035 

10 

3-0 

1.77 

I.05 

67 

F68-0036 

10 

2.3 

1.59 

1.11 

68 

F68-0037 

10 

6.2 

4.56 

1.64 

69 

F68-0038 

10 

10.3 

6.47 

I.56 

70 

F68-0039 

10 

6.9 

4.58 

1.64 

71 

F68-0040 

10 

10.5 

7.45 

1.84 

72 

f68-oo4i 

10 

10.6 

5.85 

1.71 

73 

F68-0042 

10 

15.2 

7.30 

1.42 

74 

F68-0043 

10 

15.8 

9.00 

I.76 

75 

F68-0044 

10 

15.5 

8.39 

1.4l 

7  6 

F68-0045 

10 

11.5 

7.67 

1.69 

77 

F68-0046 

12 

3-3 

6.29 

1.70 

78 

F 68-0047 

12 

3.4 

6.43 

1.30 
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oo  I  H 

l 


H  O  -4- 
r-i  cy  op 


O  UNON^J- 
CM  H  H  CO  ON 

•  •  •  •  i 

O  Q  O  O  O 


•  •  • 

O  O  Cj> 


NO  O  H 
OO  H  O 


ON  CO  CO  O 
OO  H  CM  S- 


O  0ONO 
CM  NO  OO 


9h  9  • 


QH-JHO 
jg-  jj-  j-no  c— 

•  •  •  •  • 

OOOOO 


Q  H  O 
UNOO  NO 


O  O  O 
OO  -d-  CM 


S<S  8 


CO  CM  ‘A 
Sri  O 


O  H  CM  O  CM  oo  HH00-3-  H  oo4 


ONOOOOO-d- 
CM  CM  NO  3-  NO 

•  •  •  •  • 

O  O  O  H  H 


OO  loss 
ON  S-  H  ON  H 
•  •  •  •  ( 

OO  IA  o  IA  ON 
H  H  H 


IA  0O  OO 
NO  CO  IA 


3ft& 

•  •  • 

NO  H  NO 


NO  O  ON 
NO  OO  H 


3SS8 


-3-  CO  On 
-3- CO  o 

•  •  • 

tA  O  NO 

H  H 


IA  ON -3- 
-3-  CO  CO  CM 
•  •  •  • 
•at  IA  CO  CM 


CO  ON  ON 
On  -d-  NO 


O  H  CM  O  CM  OO  H  CM  CM  OO  H  CM  CM 


NO  CO 


NONONONONO  COCOCO  OOO 


VO  VO  vo  oo  oo 
rH  rH  pH  rH  pH 


Table  C8 


Summary  of  Predicted  Data  and  Comparison  of  Prediction 

Accuracy  for  M37  Truck 

Obstacle 

Height 

in. 

Impact 

Speed 

mph 

Driver's  Seat 

Peak  Vertical  Acceleration 

Predicted 

g's 

Measured 

g's 

Deviation 

g's 

Percent 

Error 

6 

5.0 

0.70 

1.40 

-0.70 

-50 

6 

10.0 

0.85 

2.98 

-2.13 

-71 

6 

13.0 

0.94 

3.72 

-2.78 

-75 

6 

20.0 

0.85 

— 

— 

8 

5.0 

1.41 

1.64 

-0.23 

-14 

8 

7-5 

2.43 

4.10 

-1.67 

-41 

8 

10.0 

4.00 

5.70 

-1.70 

-30 

8 

11.0 

4.29 

5.92 

-1.63 

-28 

10 

2.5 

1.78 

1.60 

0.18 

11 

10 

5.0 

2.52 

2.94 

-0.142 

-14 

10 

10.0 

5.03 

7.00 

-1.97 

-28 

12 

2.5 

1.90 

5.90 

-4.00 

-68 

12 

5.0 

3.40 

— 

— 

— 

12 

10.0 

6.09 

— 

— 

— 

PLATE  Cl 


VEHUiCftL  *.rCELCWJlf  IQM  XT  Qltl'r  ETf  1 3EAT 


PLATE  C2 


PLATE  C3 


PEAK  VERTICAL  ACCELERATION  AT  DRIVER’S  SEAT,  g's 


IMPACT  SPEED,  MPH 
b.  PREDICTED 


SUMMARY  OF  MEASURED  AND  PREDICTED  SPEED- 
OBSTACLE  HEIGHT-PEAK  VERTICAL  ACCELERA- 
TION  RELATIONS,  M37  TRUCK 


PLATE  C4 


i 


IMPACT  Sf^CED, 


PREDICTED  PEAK  VERTICAL  ACCELERATION,  g  S 


PLATE  CO 


PREDICTED  SPEED,  MPH 


MEASURED  SPEED,  MPH 


NOTE:  NUMBERS  NEAR  PLOTTED  POINTS  INDI¬ 
CATE  OBSTACLE  HEIGHT  IN  INCHES. 

LEGEND 


&  M60A1  TANK 
O  M37  TRUCK 


COMPARISON  OF  MEASURED  AND  PREDICTED 
SPEEDS  AT  WHICH  PEAK  VERTICAL 
ACCELERATION  EQUALS  2.5  g’s 


PLATE  C7 


MEASURED  PEAK  LONGITUOINAL  ACCELERATION  AT  DRIVER’S  SEAT,  g's 


IMPACT  SPEED.  MPH 


0  S  10  15  20 


OBSTACLE  HEIGHT,  IN. 

LEGEND 

O  -  MEASURED,  2.0  g 
•  -  PREDICTED,  2.0  g 
X  -  MEASURED,  1.0  g 
®-  PREDICTED,  1.0  g 

SPEED-OBSTACLE  HEIGHT  RELATIONS 
FOR  M60A1  TANK  AT  1.0-fl  AND  2.0-g 
PEAK  LONGITUDINAL  ACCELERATIONS 


PLATE  CM 


COMPUTER  PROGRAM  FOR  8MJLAIIHG  DYNAMIC  RESPONSE  OF  M6QA1  TANK,  AND 
DICTIONARY  OF  PROGRAM  VARIABLES 


Pro  term 


TANK 

1SLIB# 

DIFFEQ 

2SLIB#  ALGEBR 

3STTY# 

120 

4SRPC 

5SNDM 

6SSAV 

100 

COMMON 

110 

COMMON 

120 

COMMON 

130 

COMMON 

140 

COMMON 

150 

COMMON 

160 

COMMON 

170 

COMMON 

1306 

VARI 

190 

COMMON 

200 

COMMON 

210 

COMMON 

220* 
230 
240 
250 
260 
270 
230 
29  0 
300 
310 
320 
330 
340 
350 
360 
370 
330  A 
390 
400 
41  0C- 


F3RCD1 »  F0RCD2#  F0RCD3.  FORCD4#  FORCD5#  F0RCD6 
FORCWl#  F0RCW2#  FORCW3#  FORCW4#  F3RCW5#  F0RCW6 
FORCKl#  F0RCM2#  F0RC3 3#  FORCN4#  FORCH 5*  F0RC36 
FORCrfl#  F0RCK2#  F0RCK3#  FORCK4#  F0RCK5#  F0RCK6 
FORTI 1#  FORCT1#  F3RCT2# F0RCT3#  F0RCT4# F0RCT5 
SPDEFl#  SPDEF2# SPDEF3# SPPEF4# SPDEF5# SPDEF6 
DSPDF1 #  DSPDF2# DSPDF3# DSP  DF4# DSPDF5# DSPDF6 
VAR1*  VAR2#  VAR3#  VAR 4*  VARS#  VAR6#  VAR7.  VAR8,  VAR9#  VAR10# 

MR}  • »  VAR  1 2.  VAR  1  3#  VAR  1  4#  VAR  1  5#  VAR  1  6#  VAR  1  7#  VAR  18 ,  VAR  1 9 

DR  VI  *  ZETA1  #H3RZ1«£TA1#AXL1  1#  AXL21#  AXL31#  AXL41*  AXL51#  AXL61 
DRV2#  Z  ETA2#  R0RZ2#  ETA2#  AXL  12#  AXL22#  AXL32*  AXL42#  AXL52#  AXL62 
DDRV2#  DZETA2#  D33RZ2,  DETA2#  DAXL  12#  DAXL22#  DAXL  32#  DAXL42# 
DAXL52#  DAXL62 
COMMON  3#  N STEPS#  HO RMOM 

COMMON  THRESHC  72),  GAMMAC72),  SI  GMAC  72) # SEGDEFC 72) # YC 8 6) 

COMMON  THC  4) 

DIM ENSI ON  FORCWC 6) #  FORCKC  6) #  SPDEFC  6) . DSPDF(  6) 

FORCHC  6)  #  FORCDC  6) 

DISPLC  10)#  VELCTY (10)#  ACCI  SSC  10),  ACCGS(  10) 

ACCMAXf 10)#  ACCMINC 10),  SACCSQC 1 0)  ,  RMSACCC 1 0) 

FI  DC  I2)#VARID(  10) 

<  FORCH 1#  FORCHC 1 >>#( FORCD1 # FORCDC 1 ) ) 

C  FORCW1#  FORCWC  1)),C  FORC* l#  FORCKC  D) 

C  SPDEF I #  SPDLFC 1 ) ) # ( DSPDF1# DSPDFC 1)> 

C  DRV  1#  DISPLC  1  )  )  #  C  DRV2#  VELCTY  C  P) 

C  DDRV2.  ACCI  SSC 1 ) > 

I  BELL/453  752/ 

VARI  U/5HV#  DRV#  5HV#  C-  G#  5HH#  C-  G#  5H  PITCH#  5HAXLE1#  5HAXLE2# 
5HAXLE3#  5HAXLE4#  5HAXLE5#  5HAXLE6/ 

PRI  NT#  "MURPHY  'S  M- 60  TANK  MODEL" 

PRINT# 

. . --DATA  INITIALIZATION . 


DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
EQUI  VALENCE 
EOUI  VALENCE 
EQUI  VALENCE 
EOUI  VALENCE 
EQUIVALENCE 
DATA 
DATA 


423 

CALL  3PENFC 1# "THRESH") 

430 

REA  DC  1#  )  C  THRESHC  I  ),  1=  1 

#  72) 

440 

CALL  CLOSEFC1) 

450 

CALL  OPENFC 1# "GAMMA") 

460 

REA  DC  1,  )  C  GAMMAC  I  )#  1=  1, 

72) 

470 

CALL  CL3SEFC 1) 

430 

CALL  OPENFC  1,  "SI  GMA") 

490 

REA  DC  l#)CSIGMACI),Isp 

72) 

500 

CALL  CLOSEFC1) 

510 

THC 1)= 13. 

520 

THC  2> -  1  0. 

530 

THC  3)  *  8. 

(1  of  15  sheets) 


TAM*  COMTIMUEO 


540  T*C  4)  »  3. 

550  DO  10  1*1.6 

560  FORCWC l )  a  0. 

570  FORC*< I )  =  0. 

530  SPDEFC  I  ) a  0. 

590  10  DSPDrCI)a0. 

600  DO  20  1*1#  10 

610  DISPL(I>=0. 

620  VELCTYCI>*0. 

630  ACCI SSC I ) *  0« 

640  ACCGSCI>*0. 

650  RMSACC(I)*0. 

660  ACCM AXC I )  =  0. 

670  ACC* I. MCI  )«0. 

630  20  SACCSOC  Ds0i 

690  DO  30  1*1*36 

700  30  YU )=0. 

710  ZETA1  a  -5.79 

720  *ORZ1=0. 

730  ETA  1  a  -.0039 

740  AXL  0*  1  3  • 

750  AXL  11  a  -.966 

760  AXL21  *  -.970 

770  AXL31  =  -.942 

730  A  XL  41  =  -.9  13 

790  AXL51  *  -.334 

300  AXL61  =  -.356 

310  TP=  0 • 

320  TIP*.  5 

330  T=0. 

340  DELTAL*33./14. 

3  50  MPL*  4 

360  MSTOP=0 

3  70  JJa  1 

330  **.001 

39 0C . DATA  READ  IN- . 

900  PRINT."GIVE  THE  FOLLOWING  I  NFORM  ATI  ON" 
910  PRINT. 

920  PRINT.  "NAME  OF  PROFILE  INPUT  FILE" 

930  READ  1.  FI  NAME 

940  PRI  NT.  "TAN*  VELOCITY  IN  M.P.*." 

9  50  READ.  XM PI 

9  60  PRINT.  "PRINTOUT  TIME  INTERVAL" 

9  70  READ.  TIP 

930  PRI  NT.  "NAME  OF  OUTPUT  DATA  FILE" 

990  READ  1.FNAME2 

1000  PRINT.  "*****  EN  D  OF  INPUT  DATA*****" 

1010  VEL*XMP**  17«6 

1020  DELTAT*  DELTAL/VEL 

1030  NSTEPS* DELTAT/* 


(2  of  15  sheets) 


TANK  CONTI  CUED 


1040 
1050 
1060 
1070 
1080 
1090 
1  100 
1110 
1 120 
1130 
1 140 
1150 
1  160 
1170 
1  180 
1190 
1 200 
1210 
1220 
1230 
1240 


H=DELTAT/NSTEPS 

PRINT  2  *  XM  Prl  *  VEL*  DELTAL*  DEL  TAT*  NSTEPS*  8 
WRI TEC  2J  1  1 ) 

CALL  OPENFC  1*  FI  NAME) 

READ!  1*  DFID 
PRINT* 

PRINT*  "INPUT  PROFILE  IS«" 

PRINT  DFID 
PRINT  3 

WRI  TEC  21  12)  FID*  FI  NAME 

WRI  TEC  2*  13>XMPr|*  VEL*  DELTAL*  DEL  TAT*  NSTEPS*  M 
WRI  TEC  21  14)  VARI  D 
GO  TO  19  0 

150  IFCJJ-2)  50*  40*  50 
50  READC 1* 15)YINPUT 
CALL  EOFTSTC 1* JJ) 

GO  TO  C  60*  40)  JJ 
40  NSTOP*NSTOP+ 1 
60  1*85 
70  YCI+lDrci) 

1*1-1 


1250 

1260 

1270 

1230 

129  0 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 


I  FC  I  >70*80*70 
80  YC  1)=T INPUT 
CALL  DIFFEQ 
NPL*NPL+ 1 
T“  T+DELTAT 
DO  90  1*1,  10 

90  ACCGSCI)*ACCISSCI>/336. 

ACCGSC  4) * ACCI  SSC  4) 

DO  100  I*  D  10 

ACCMAXCI)*AMAX1C  ACCMAXCI  ),  ACCGSCI  )) 
ACCMINCI )* AMIN  1C  ACCMINCI ),  ACCGSCI) ) 
SACCSQC  I  )*SACCSOC  I  >«■  ACCGSC  I  )  *ACCGSC  I ) 
1 00  RM  SACCC I ) *  SQRTC  C SACCSQC I ) *  DEL  TAT> / T) 
190  IFCT-TP)  1  10,  120,  120 


139  0  120  TP=TP+TIP 


1400 
141  0< 


PRINT  4*  T»  YC  1)  *  C  VARI  DC  I  ),  DI  SPLC  I  ), 
RMSACCCI)*I=D  10) 


VELCTYCI ) *  ACCGSC  I ) 


1420  110  WRITEC2*  16)T  YC  1DDISPL*  VELCTY 
1430  IFCNPL- 1 1)  140>  130*  140 
1440  130  WRITEC2* 1^)VARID 


*  ACCGS* RMSACC 


1450  NPL*0 


I 


140  IF  C  NSTOP-36)  150*  160.  150 

160  *5><  VARID(I)*  ACCMAXCI ),  ACCMINCI),  I  *  d  10) 

PRINT  5* C VARI  DC  I  )* ACCMAXCI), ACCMINCI), I* D 10) 

I  FC  FNAME2-  6HN0FIL  E)  170*  19U*  1  70 
170  CALL  CLOSEFC 2,  FNAME2,  2) 

PRINT  6*  FNAME2 


1520  180  PRINT  7,  Cl  BELL*  1*1,  40) 
1530  CALL  EXIT 


(3  of  15  sheets) 


TANK  CONTI  SUED 


1  540  1  FORMAT?  12A6) 

1550  8  FORMAT?  ///» “VELOCI  TY-“*  F5*  2*  "  MPK  ?“»F6.1»“  IPS  )“»/» 

1560*  “DELTA-L*"*  F5.  3*  3X.  “DELTA-  T«“»  F6.  4*  /» 

1570*  "NSTEPS““»  I  At  A X*  “H«“»  F7. 6) 

1580  3  FORMAT? /////2X»  4HTIME*  3X»  4HY?  I) *  12X.  5HDISPL*  5X»  5HVEL0C. 

1590*  5X»  5HACCEL*  4X»  6HRMSACC*  /) 

1600  4  FORMAT? F8.4»F7. 2*  2X»  AS*  4G10.3*  /»9?  1  7X#  AS.  4G10.3* /)) 

1610  5  FORMAT?  3 7H PEAK  ACCELERATIONS  MAXIMUM  MINIMUM*  /» 

1620*  10?  10X.A5.2X.  2FI0.3*/)) 

1630  6  FORMAT?  38HCALL  2530  TO  OBTAIN  A  LISTING  OF  FILTI.X.A6) 

1640  7  FORMAT? 40A1) 

1650  11  FORMAT?//*  37X*  45?  1H*)  */.  37X.  1H**43X.  1H*»/»37X» 

1660*  45H*  MURPHY  'S  M-  60  TANK  PROGRAM  OUTPUT  FILE  *»/» 

1670*  37X* 1H*» 43X* 1H*» /. 37X* 45? 1H*)»//) 

1680  12  FORMAT? 16HIN PUT  PROFILE  I  S* X.  1 2A6*  X.  1  HI  FILE  NAME  .A6.1HJ*//) 

1690  13  F0RMAT?9H VELOCI  TY-FS.  2*  1  7H  MILES  PER  HOUR  ?F6.1» 

1700*  19H  I  NCHES  PER  SECOND)  4X*  8HDELTA-L-F5.  3*  7H  I  NCHES4X* 

1710*  8H DELTA- T-F10.8.8H  SECONDS//* 

1720*  35HNUMBER  OF  STEPS  IN  RK  G  INTEGRATION-1  4*  4X» 

1730*  12HSTEP  SIZE  H-F10.8) 

1740  14  FORMAT?//. 2X»  4HTIME3X.  4HY?  1)  12X*  10?A5»  3X)»//) 

1750  15  FORMAT?  E20.  10) 

1760  16  FORMAT?  F8  »  4*  F5«  1 »  2X»  6HDI  SPL  •  2X*  1 0F8  •  3*  /*  I  5X»  8H  VELOCI  TY » 

1770*  1  0F8  •  3*  /*  15X.  6HACCEL.2X*  10F8.3*/*  1  5X*  8HRMS.  ACC.  1  0F8  •  3* /) 

1780  END 

2000  SUBROUTINE  DIFFEQ 
2010  RH-l./H 
2020  INDEX-0. 

2030  1  INDEX-INDEX* 1 
2040  VAR1-ZETA 1 
2050  VAR2-ZETA2 
2060  VAR3-H0RZ 1 
2070  VAR4-HORZ2 
2080  VAR5-  ETA  1 
2090  VAR6-  ETA2 
2100  VAR19  -  AXL0 

2110  VAR7- AXL 1 1 
2120  VARS- AXL 12 
2130  VAR9- AXL21 
2140  VAR10-AXL22 
2150  VAR1 1-  AXL 31 
2160  VAR 12- AXL 32 
2170  VAR 1 3= AXL 41 
2180  VAR1 4-AXL42 
2190  VAR1 5-AXL51 
2200  VAR 16- AXL 52 
2210  VAR17-AXL61 
2220  VAR18- AXL62 
2230  PZETA2-ZETA2 
2240  PH0RZ2-H0RZ2 
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TANK  CONTINUED 


2250  PETA2*  ETA2 
2260  PAXL 12* AXL 12 

2270  PAXL 22s AXL22 
2230  PAXL  32*  AXL  32 

229  0  PAXL  42s  AXL  42 

2300  PAXL 52* AXL 52 

2310  PAXL 62* AXL 62 

2320  CALL  ALGEBR 

2330  CALL  EQNS<  FK  12#  FK  1  4#  FK  1 6*  FK  13*  FK  1 1 0*  FK  1  1 2.  FK  1 1  4*  FK  1 1  6*  FK  1 13) 

2340  FK 1 1*H* VAR2 

23S0  FK 13*H*VAR4 

2360  FK  1 5*rl*  VAR6 

2370  FK 1 7*rl*  VAR3 

2330  FKl9*H*VAR10 

2390  FK  1  1  l»rl*VARl  2 

2400  FK 1 1 3*H*VAR1 4 

2410  FK 1 1 5*rl*VARl  6 

2420  FK  1 1  7*rl*VAR18 

2430  VAR1*ZETA1  +  FK1 1*.  5 

2440  VAR2*ZETA2+FK12*.5 

2450  VAR3*H0RZ 1 +  FK  1  3*  •  5 

2460  VAR4*H0RZ2+FK14*.5 

2470  VAR5*ETA1  +  FK15+.  5 

2430  VAR6*  ETA2+  FK 1 6* • 5 

249  0  VAR7* AXL  1  1  ♦  FK  1  7*  •  5 

2500  VAR8- AXL 1 2+ FK 1 3  * • 5 

2510  VAR9*AXL21+FK19*. 5 

2520  VAR 1 0*  AXL22+  FK 1 1  0*  • 5 

2530  VAR1 1*AXL31+FK 1 1  !*•  5 

2540  VAR1 2* AXL  32* FK 1 1 £*  •  5 

2550  VAR1 3* AXL 41* FK 1 1 3* • 5 

2560  VAR 1 4*  AXL  42*  FK 1 1  4*  • 5 

2570  VAR1  5*AXL51  +  FK  1 1  5*  •  5 

2580  VAR1 6* AXL  52+FK 1 1 6* • 5 

2590  VAR1  7*AXL61  +  FK  1  1  7*  •  5 

2600  VAR 13* AXL 62+ FK1 13*. 5 

2610  CALL  ALGEBR 

2620  CALL  EON  S<  FK 22*  FK2  4*  FK26*  FK 23*  FK  2 1  0*  FK 2 1  2*  FK 2 1  4*  FK  2 1  6*  FK  2 1 3  > 

2630  FK2 1*H+ VAR2 

2640  FK23=rl+ VAR4 

2650  FK25*H*VAR6 

2660  FK27*H* VAR3 

2670  FK29*H* VAR1 0 

2630  FK21  l*rl  +  VAR12 

2690  FK21  3*rC  +  VARl  4 

2700  FK21 5*K*VAR1 6 

2710  FK21 7*H* VAR13 

2720  VAR1»ZETA1+.29239322*FK21+.20710673*FK1 1 

2730  VAR2*ZETA2+ • 29  239  322*FK22+ • 2  0  7  1  06  73*  FK 1 2 
2740  VAR3*H0 RZ 1+. 29239 322*FK23+. 2071 0673* FK 13 
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TAM*  CONTINUED 


27S0  VAR4=H0RZ2*.29289322*FK24*  <20710678*FK 14 

2  760  VARS*  ETA 1 * . 29 289  322*  FK 25* . 207 1 06  78*  FK 1 5 

2770  VAR6*  ETA2* • 29 289322* FK26*«  207 1 0678* FK 1 6 

2780  VAR7*AXL  1 1* • 29289 322*FK27*«  2071 0678* FK 1 7 

2790  VAR8*AXL  12* . 29 289 322* FK23* . 20710678* FK 18 

2800  VAR9*AXL2 1* • 29  28  9  322*FK29* • 2071 0678* FK 19 

2810  VAR  1 0* AXL 22* • 29  239  322* FK 21 0*  *207 1 0673* FK 1 1 0 

2820  VAR1 1*  AXL31* • 29  289  322* FK2 1 1*  »207 1 0678*  FK 1 1 1 

2830  VAR 12* AXL 32* • 29  239  322*  FK21 2* • 207 1 0678*FK 1 12 

2840  VAR1 3* AXL4 1* • 29289  322*  FK21 3*«  207 1 0678*FK 1 1 3 

2850  VAR1 4= AXL 42*. 29 289 322* FK 21 4*.20710678*FX 1 1 4 

2860  VAR1 5* AXL  5 1* • 29  239  322*FK2l 5* . 2 071 0678* FK 1 1 5 

2870  VAR 16* AXL 52* • 29289  322*FK21 6* *207106 78*FK1 16 

2330  VAR  1 7* AXL 6 1* • 29239  322* FK 2 1 7* • 2071 0678*  FK 1 1 7 
289  0  VAR 1 3*  AXL  62* • 29  289  322*  FK218*.20710678*FK118 

29  00  CALL  ALGEBR 

2910  CALL  EQN  S C  FK  32*  FK  3  4#  FK  36#  FK  33 #  FK  3 1  0*  FK  3  1 2 j  FK  3 1  4*  FK  3 1  6#  FK318) 

2920  FK31*K*VAR2 

2930  FK33*K* VAR4 

29  40  FK35=K*  VAR6 

29  50  FK37*K*VAR8 

2960  FK39*K* VAR  1 0 

2970  FK3U  =  K*VAR12 

2980  FK31 3*K* VAR1 4 

2990  FK31 5*r|*VAR16 

3000  FK31 7* H* VAR  18 

3010  VAR1*ZETA1- . 7071 0678* FK2 1* 1 . 7071 0678* FK 31 

3020  VAR2*ZETA2- . 70710678*FK22* 1 . 7071 0678* FK 32 
3030  VAR3*HORZ 1- . 70710678*FK23* 1 . 7071 0678* FK33 
3040  VAR4*K0RZ2-  .  7  0  71  06  78*FK24*  1 . 7071  0678*FK34 

3050  VAR 5* ETA  1- • 707 106 78* FK 25* 1 • 7 071 0678* FK 35 

3060  VAR6*ETA2- . 70710678*FK26* 1 . 7071 0678*FK36 
3070  VAR7*AXL1 1-.  7071  0678*FK27*1. 70710673*FK37 
3080  VAR3*AXL 12- • 7071 0673* FK28* 1 . 7071 0678*FK38 
309  0  VAR9* AXL2 1“ • 70  71 0678* FK29* 1 . 7071 0678* FK  39 
3100  VAR10=AXL22- • 707 1 0678* FK2 1 0* l . 7071 0673*FK3 1 0 

3110  VAR1 1*AXL31-«7071 06 78* FK 2 l 1+ 1 . 7071 0673*FK31 1 

3120  VAR  12* AXL 32- • 707 1 0678*FK2  1  2* 1 . 7  0  71  06  78* FK 312 

3130  VAR13=AXL41-. 7071 0678* FK2 l 3* l . 7071 0678* FK 31 3 

3140  VAR1 4= AXL 42“ • 7071 0678* FK 2 l 4* 1 . 7071 0678*FK31 4 
3150  VAR 1 5* AXL 51- . 707 1 0673*FK 2 1 5* l . 70710678*FK3l 5 
3160  VAR1 6* AXL  52- • 7071 0673* FK 2 1 6* 1 • 7071 0678* FK  3 1 6 
3170  VAR1 7=AXL61-  •  7071 0678*FK2l 7* 1 • 70  71 06  78* FK 31 7 
3130  VAR18*AXL62-. 707 1 0678* FK 2 13* 1. 7071 0673+FK318 

319  0  CALL  ALGEBR 

3200  CALL  EQNS(  FK  42*  FK44*  FK  46#  FK  43»  FK  41  0»  FK41  2»  FK41 4»  FK  41  6»  FK  418  > 

3210  FK4l  =  K*  VAf<2 

3220  FK43*K* VAR4 

3230  FK45=H*VAR6 

3240  FK47*K* VAR8 
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TANK  CONTINUED 


32S0 

3260 

3270 

3280 

329  0 

3300 

33101 

3320 

33301 

3340 

33S01 

3360 

33701 

3380 

33901 

3400 

34101 

3420 

34301 

3440 

34501 

3460 

34701 

3430 

34901 

3500 

35101 

3520 

35301 

3540 

35501 

3560 

35701 

3580 

359  01 

3600 

36101 

3620 

36301 

3640 

36501 

3660 

3670 

3630 

369  0 

3700 

3710 

3720 

3730 

3740 


FK49  =  H*VAR10 
FK41 1«H*VA»12 
FK  41 3=rf*VAR 1 4 
FK415*H*VAR16 
FK41 7=H* VAR18 

ZETA 1=ZETA 1* •  1 66667*FK 1 1*. 09  763107* FK 2 1* . 569  0  3  5  59* FK 31 
+  •  1 66667*FK  4 1 

ZETA2=Z  ETA2* •  1 66667*FK 12*. 09  7  6  31  07* FK 22+ • 569 03 559 *FK 32 

♦  .  166667*FK42 

H3RZ l  =  HORZl*.  1 66667*FK 13*. 09  7631 07*FK23+ • 5 69 03 559 *FK 33 
166667*FK43 

H0RZ2=.H0RZ2*  <  1  66667*FK  1  4*  <  09  7631  07*FK24*  .  569  03559*FK34 

♦  • 1 66667*FK  44 

ETA  1  =  ETA  1*.  1 66667*FK 1 5* • 09  7  6  3  1  07* FK 2 5+ • 569  03559* FK 35 
♦. 166667*FK45 

ETA2=ETA2*.  1 66667* FK  1 6* •  09  7631  07*FK2  6+.  569  03559*FK36 
♦. 166667*FK46 

AXL1  1sAXl.11*.  1  666 67* FK  1  7*  •  09  7631  07*  FK  27*  •  56903559*  FK  3 7 
166667*FK47 

AXL  1  2=  AXL  1  2*  •  1  66667* FK  18*.  09  7  6  31  07*FK28*  .  56903559*FK38 
*•  1 66667* FK  48 

A.\L  2  1  =  AXL  2  I  *  .  1  6666  7*  FK  19*.  09  763 1 07*  FK  29*.  56903559*FK  39 
*•  166667* FK  49 

AXL22=AXL22+. 166667*FK1 1 0+ . 09 7631 07* FK21 0* . 569 03559*FK 31 0 
1 66667*FK  4 1 0 

AXL  31  =  AXL  31*.  16t,  .67*FK1  1  1*  .  09  7  6  3  1  07*FK2 1  1  *  .  569  0  3559*  FK  3 1  1 
* • 166667* FK41 1 

AXL 32* AXL 32* . 1 66667* FK 1 1 2* • 09  7  6  3  1  07*  FK21 2* • 569  0  3  559*FK31 2 
*. 166667*FK412 

AXL 41  =  AXL 41*.  1 6666 7*FK 1 1 3* • 09 7631 07* FK2 1  3* . 569 03559* FK 31 3 
16666 7* FK 413 

AXL 42=  AXL 42*.  1  6666 7*FK  1  1  4*  •  09 7631  07* FK21  4*  •  569  0  3  5  59*FK3  1  4 
+• 166667*FK414 

AXL 51  =  AXL 51*. 1 66667*FK 1 1 5*. 09  7631 07* FK2  1  5* . 569  03559 *FK 31 5 
♦. 166667* FK415 

AXL  52* AXL  52* •  1 66667*FK 1  1 6* . 09  76  31  07* FK 2  1  6* . 569  0  3  559*FK  31 6 

♦  •  1 66667*FK  4 1 6 

AXL  61  =  AXL  6 1*  •  1  6666 7* FK  1 1  7*  .  09  763 1  0  7*FK2 1  7*  •  569  03559 *FK 31  7 
♦. 166667*FK417 

AXL 62= AXL  62* .  166667*FK 1 18*. 09  7631 0  7*FK218* . 56903559* FK 31 8 

♦  •  166667*  FK413 

DZ  ETA2*  <  Z  ETA2- PZ  ETA2 ) *  RH 
DH0RZ2=  <H0RZ2-  PHORZ2)  *RH 
DETA2=< ETA2- PETA2) *RH 
DAXL  1 2=  C  AXL  1 2-  PAXL  1 2)  *  RH 
DAXL  22=  <  AXL 22-  PAXL  22)  *  RK 
DAXL  32=  <  AXL  32-  PAXL  32)  *  RH 
DAXL  42=  <  AXL  42-  PAXL  42  )  *  RH 
DAXL  52=  <  AXL  52-  PAXL  52)  *  RH 
DAXL  62=  <  AXL  62- PAXL  62) *RH 
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TANK  CONTI  NUED 


3750  DRV1  =  ZETA1  ♦  25  *  ETA  1 

3760  DRV2  =  ZETA2  ♦  2S  *  ETA2 

3770  DDRV2  =  0ZETA2  ♦  25  *  DETA2 

3780  I  F!  I  NDEX-NSTEPS) 1 »  2#  2 

3790  2  RETURN 
3800  END 

3310  SUBROUTINE  EQNSC  FK2#  FK  4#  FK6#  FK8#  FK  I  0,  FK  1 2#  FK  I  4#  FK  I  6*  FK  13) 

3320C 

3830C . VERTICLE  C- G  EQUATION . 

3840  FK  2= H  *  I  C  -  FO  RCK  1  -  FO  R  CK  2-  FO  RCK  3-  F  0  R  CK  4-  FO  R  CK  5-  FO  RCK  6 
33S04  - FORCDI- F0RCD2" F0RCD3” F0RCD4- F0RCD5- F3RCD6) *«  008“ 38  6.) 

33  60C 

33  70C . HORIZONTAL  C- G  EQUATION . 

3830  FK4=H*I  FORCH  1  ♦F0RCH2*  FORCn  3* FORCH  4*  FORCH  5*  FORCH  6)  * .  008 
339  0C 

39  0  0C . MOMENT  AT  C-  G  EQUATION . 

39 10  FK6=H*I- 77. +FORCK1- 44.* FORCK 2- 1 1  .  *FO RCK 3+22. *FORCK 4 
39  204  ♦ 5S«  * FORCK  S+ 83. *FO RCK 6- 77.* FORCDI- 44.*  F0RCD2- 1 1 .*F3RCD3 

39  304  ♦  22.*F0RCD4*S5.*F0RCD5'*'33  •  *F0RCD6-H0RM0M>  /58  1  700. 

39  40C 

39  50C . BOGIE  EQUATIONS  C6> . 

3960  FK3=H*I FORCK 1+FORCD1* FORCT1- FORT 1 1+FORCW1- 1 420. >*.271 7 
39  70  FK 10=H*I F0RCK2+F0RCD2- FORCT1 *F0RCT2* F0RCW2-  1420.) *.271 7 
398  0  FK 1 2«H*I  FORCK  3*F0RCD3- F0RCT2*F0RCT3*F0RCW3- 1 420. ) *. 271 7 
399 0  FK 1 4=H*I FORCK 4* F0RCD4- F0RCT3*F0RCT4*F0RCW4-  1 420. >  * . 27 1 7 
4000  FK  1  6=  H  *  I  F0RCK5*F3RCDS“F3RCT4*F0RCT5*F0RCW5-  1420.)  *.271  7 
4010  FK  18*H*I  FORCK  6*  F0RCD6- F0RCTS*F',RCW6“  1420.  >*.271  7 
4020  RETURN 
4030  END 

5000  SUBROUTINE  ALGEBR 
5010  DO  95  I  *  1*  6 

5020  FORCWII)  =  0. 

5030  95  FORCH C I )  =  0. 

5040  VAR  19=  0. 

5050  DO  20  1= 1*  4 

5060  I  FIYII  >-THII  >>20#  20#  30 

S070  30  VTH= YI I > - THC I ) 

5080  VAR19=AMAX1( VAR19# VTH) 

509  0  20  CONTINUE 

5100  DO  100  Is  1#  12 

5110  SEGDEFII )  *Y  1 1  ♦ 4) -  VAR 7- THRESH  1 1  ) 

5120  I  FC SEGDEFII ) )  50# 60# 60 
5130  50  SEGDEFII  )  =  0. 

5140  60  FORCWl=FORCWl* SEGDEFII  )* GAMMA! I  ) 

5150  100  FORCH  1  =  FORCH  1  ♦  SEGDEFt  I )  *SI  GMAII) 

5160  DO  200  I  =  1  3#  2  4 

5170  SEGDEFII  >  =  Y<  I  ♦ 6) - VAR9- THRESH!  I  ) 

5180  IFISEGDEFII))  150#  160#  160 
5190  150  SEGDEFI I )  =  0. 

5200  160  F0RCW2=F0RCW2* SEGDEFII )*GAMMAII) 
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tank  continued 


5210  200 
5220 
5230 
5240 

5250  250 
5260  260 
5270  300 
5280 
529  0 
5300 

5310  350 
5320  360 
5330  400 
5340 
5350 
5360 

5370  450 

538  0  460 

539  0  500 
5400 
5410 
5420 

5430  550 

5440  560 

5450  600 

5460 

5470 

548  0 

5490 

5500 

5510 

5520 

5530 

5540 

5550 

5560 

5570 

5580 

5590 

5600  720 
5610 

5620  710 
5630  730 
5640 

5650  740 
5660  700 
5670 
5680 

569  0  810 
5700 


FORCH 2= FORCH 2*  SEGDEFCI)*SI  GMA(  I ) 
DO  300  Is 25*  36 
SEGDEFC I)=YCI*8>-VARl 1-THRESHC I ) 

I  FC  SEGDEFC  I ) )  250#  260#  260 
SEGDEF( I )=  0. 

F0RCW3=F0RCW3*  SEGDEFC I ) *  GAMMA!  I ) 
FORCH 3* FO RCH 3* SEGDEFC I  )*SI  (MACD 
DO  400  I  =  37#  48 

SEGDEFC  I  >  =  Y(I*  1  0)  -  VAR  1  3-THRESH  C  I  > 
I F( SEGDEFC I ) ) 350#  360#  360 
SEGDEFC I )s0» 

F0RCW4*  FORCW4+  SEGDEFC  1  ) *  GAMMAC  I  ) 
FORCH 4* FORCH 4*  SEGDEFC  I  )  *  SI  (MACI  ) 
DO  500  I  =  49.  60 

SEGDEFC I )aYCI* 1  2) -VAR  1  5- THRESH  C  I > 
I FC  SEGDEFC I ) ) 450#  460#  460 
SEGDEFC I ) s  0. 

F0RCW5=  F0RCW5*  SEGDEFC  I  )*GAMMAC  I  ) 
FORCri 5*  FORCH  5*  SEGDEFC  I  )  *SI  GMACI  ) 
DO  600  I*A|.»  72 

SEGDEFC!  ) aYC  !♦  14)-  VAR  1  7- THRESH  Cl  ) 
I FC  SEGDEFC  I  ) ) 550#  560#  560 
SEGDEFC I )=0. 

F0RCW6aF0RCW6*  SEGDEFC I ) ♦ GAMMAC  I  ) 
FO RCH 6* FORCH  6*  SEGDEFC  I  )  *  SI  Ol  AC  I  ) 
SPDEF1*  VAR1+  77«*  VARS- VAR 7 
DS  PDF 1 s  VAR2+ 7  7.*  VAR6- VAR8 
SP0EF2* VAR 1*44.  * VAR5-VAR9 
DSPDF2*  VAR2+ 44. *  VAR6- VAR 1 0 
5PDEF3* VAR1* 1 1 .*VAR5-VARl  1 
DSPDF3*  VAR 2* 11.*  VAR6- VAR1 2 
SPDEF  4s  V AR 1 - 22 . *  VAR5- VAR 1 3 
DSPDF4* VAR2- 22-*VAR6- VAR1  4 
SPDEF5* VAR1- 55. * VAR5- VAR1  5 
DSPDF5= VAR2- 55. * VAR6- VAR  1  6 
SPDEF6s VAR  1-83.* VARS- VAR  17 
DSPDF6* VAR2-88 .*VAR6-VAR13 
DO  700  I*l#6 

I  FC  SPDEFC  I  )-  .  402)  710#  71  0*  720 
SPDEFC I ) s . 402 
DSPDFC I )s0. 

I FC  SPDEFC I >♦ 12. )  730#  740#  740 
FORCKC I)= 29993 .* SPDEFC I ) *3399 72. 

GO  TO  700 

FORCKC I )» 1667. * SPDEFC I ) 

CONTINUE 
DO  800  I*l#6 

I FC  ABSC  DSPDFC I )) -  1 .) 8 1 0#  320#  3  20 
FORCDC I )»  2750.* DSPDFC  I ) 

GO  TO  8  00 
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TANK 


5710 

5720 

5730 

5740 

5750 

5760 

5770 

578  0 

5790 

58  00 

5810 

58  20 

5830 

5840 

58  50 

58  60 

58  70 

588  0 

5890 

5900 


820  FORCD( I  >»  SI GN<  2750./ DSPDF< I >  > 
800  CONTINUE 

I  F<  VAR7-  VAR  19 >9  50/  9  50/960 
950  FORTl 1* 600.*( VAR7- VAR19) 

GO  TO  9  70 
960  FORTH 


<  VAR9- VAR 7) 

(  VAR1 1- VAR9) 

( VAR 13- VAR 1 | > 
(VAR15-VAR13) 
(VAR  17- VAR 15) 


9  70  FORCT1  ■  375.  * 

F0RCT2  ■  375.  * 

F0RCT3  «  375.  * 

.  FORCT4  «  375.  * 

F0RCT5  ■  375.  * 

HORMOMa  0. 

HORMOM«FORCH 1*< 46.+SPDEF1+  77. *VAR5) 
HORMOM»HORMOM+FORCH2*(  46.  +  SPDEF2+44.  * VAR5) 

HORMOM-rfORMOM*FORCH3*(46.-SPDEF^n  IS 
HORMOM.rtORMOM*FORCK4*«46.*SPDEF4-88.IvAR5> 

H0RM0M»H0R'i0M+F0RCrf5*(  46.+SPDEF5-  55*I\/AR5) 

RESjSJ*rf0RM0M+F0RCH6*<46,  +  SPDEF,6‘88**VAR5> 

END 
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THRESH 


100  7. 4, 5. ,3., 1 .8, .6, .05, 
110  .05, .6, 1 .8, 3 ,  ,5 .,7,4, 
120  7. 4, 5. ,3., 1 .8, .6, .05, 
130  .05, .6, 1.8,3. ,5 .,7.4, 
140  7. 4, 5. ,3. ,1.8, .6, .05, 
150  .05, .6, 1 .8,3. ,5  ., 7.4, 
160  7. 4, 5. ,3. ,1.8, .6, .05, 
170  .05* .6, 1 .R,3.,5 .,7.4, 
180  7. 4, 5. ,3.,  1  .8,  .6,  .05, 
190  .05, .6, 1.8, 3. ,5. ,7. 4, 
200  7. 4, 5. ,3., 1 .8, .6, .05, 
210  .05, .6, 1.8, 3. ,5. ,7. 4 


r,4  iia 


100  2296, 2828, 3276, 3624, 3864, 390*, 
110  3984,3864,3624,3276,2828,2296, 
120  2296,2828,3276,362/1,3864,3984, 
130  3984,3864,3624,3276,2828,2296, 
140  2296,2828,3276,362/1,3864,3984, 
150  3984,3864,3624,3276,2828,2296, 
160  2296,2828,3276,3624,3864,3984, 
1 70  3984,3864,3624,3276,2828,2296, 
180  2296,2828,3276,3624,3864,3984, 
190  3984,3864,3624,3276,2828,2296, 
200  2296,2828,3276,3624,3864,3984, 
210  3984, 3864, 3624, 3276, 282«, 2296 


SIGMA 


1  00  71  1? ,6150,4997 ,3687 ,2250,757 , 

1  10  -757,-2250,-3697,-4997,-6150,-71 19, 
120  71  19, 6 150, 4997, 3 697, 2250 ,75 7, 

130  -757, -2250, -36B7, -4997, -6150, -71  10, 
140  71 19,6150,4997,3697,2250, 75 T, 

150  -757,-2250,-3697,-4097,-6150,-7119, 
160  71 19 , 6150, 4097,3697, 2250 ,757, 

170  -757,-2250,-3697,-4097,-6150,-7119, 
190  71 19, 6150,4997, 3697, 2250, 757, 

100  -757,-2250,-3697,-4097,-6150,-71 19, 
200  7110,6150,4997,3697,2250,757, 

210  -757,-2250,-3697,-4997,-6150,-71 19 


FIMAKE 


10C  THIS  PROGRAM  IS  USED  TO  MAKE  OBSTACLE  DAYA  FILES  TO  BE  INPUT  TO 
11C  THE  VEHICLE  DYNAMICS  PROGRAM. 

100  DIMENSIONXd000).FILEIDd2) 

110  PRINT,  FILE  NAME 

120  READ  1 , FI  NAME 

130  PRINT, "FILE  IDENTIFICATION 

140  READ  1 , FILEID 

1 53  PRINT, "NO.  OF  POINTS 

1 60  READ.  N 

170  PRINT2.N 

190  READ, (Xd ) , Is  1 , N) 

190  VRITEd.DFILEID 
200  WRITEd  ,3)  (X(  I  ),  1=  1  ,N) 

210  PRINT, "TYPE  NO.  OF  TRAILING  ZEROES 
220  READ, M 
230  v=0 

240  WRITEd, 3HY, 1=1, M) 

250  Nr  fH-M 

260  PRI NT4, N, FI  NAME 

270  CALLCLOSEFd  ,FINAME,2) 

290  CALLEXIT 
290  1  FORMAT ( 1 2A6 ) 

300  2F0RMAT(5HTYPE,I4, 12H  DATA  POINTS) 

310  3FORMAT(E20.10) 

320  4F0RMAT("THERE  ARE  ,,I3,  POINTS  IN  FILE  ,A6) 

330  END 
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Dictionary  of  Variables 


Variable 

FORCK(I) 

FORCW(I) 

FORCH(I) 

FORCD(I) 

SPDEF(I) 

DSPDF(I) 

FORCT(I) 

FORTH 

VAR1-VAR19 

DRV1 

DRV2 

DDRV2 

ZETA1 

ZETA2 

DZETA2 

HORZ1 

IIORZ2 

DHORZ2 

ETA1 

ETA2 

DETA2 

AXL(I) 1 

AXL(I)2 

DAXL(I)2 

H 

RH 

NSTEPS 
liORMOM 
THRESH (I) 
GAMMA(I) 
SIGMA (I) 
SEGDEF(I) 


_ Description 

FORCE  OF  ITII  SUSPENSION  SPRING 

RESULTANT  VERTICAL  FORCE  OF  SPRING  SEGMENTS  OF  ITH  BOGIE 

RESULTANT  HORIZONTAL  FORCE  OF  SPRING  SEGMENTS  OF  ITII  BOGIE 

FORCE  OF  ITII  SUSPENSION  DAMPER 

DEFLECTION  OF  ITII  SUSPENSION  SPRING 

VELOCITY  OF  ITH  BOGIE  DAMPER  (DERIVATIVE  OF  SPDEF(I)) 

TRACK  TENSION  FORCE  BETWEEN  BOGIE(l)  AND  BOGIE(l+l) 

TRACK  TENSION  FORCE  OF  FEELER  SPRING 

PAST  VALUES  OF  ALL  VARIABLES  IN  RUNGE  KUTTA  SOLUTION 

VERTICAL  DISPLACEMENT  OF  DRIVER 

VERTICAL  VELOCITY  OF  DRIVER 

VERTICAL  ACCELERATION  OF  DRIVER 

VERTICAL  DISPLACEMENT  OF  CG 

VERTICAL  VELOCITY  OF  CG 

VERTICAL  ACCELERATION  OF  CG 

HORIZONTAL  DISPLACEMENT  OF  CG 

HORIZONTAL  VELOCITY  OF  CG 

HORIZONTAL  ACCELERATION  OF  CG 

PITCH  DISPLACEMENT  ABOUT  CG 

PITCH  VELOCITY  ABOUT  CG 

PITCH  ACCELERATION  ABOUT  CG 

VERTICAL  DISPLACEMENT  OF  ITII  BOGIE 

VERTICAL  VELOCITY  OF  ITII  BOGIE 

VERTICAL  ACCELERATION  OF  ITH  BOGIE 

RUNGE  KUTTA  TIME  STF.r 

RECIPROCAL  OF  H 

NUMBER  OF  STEPS  IN  RUNGE  KUTTA  SOLUTION 
MOMENT  OF  HORIZONTAL  FORCES  ABOUT  CG 
THRESHOLD  HEIGHTS  OF  ITH  SPRING  SEGMENT 
kyCOS  i 
khSIN  i 

DEFLECTION  OF  SEGMENT  SPRING(I) 
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Variable 


Description 


Y(I) 

TH(I) 

DISPL(I) 

VELCTY(I) 

ACCISS(I) 

ACCGS(I) 

ACCMAX(I) 

ACCMIN(I) 

saccsq(I) 

RMSACC(I) 

VARID(I) 

FID 

IBELL 

TP 

TIP 

T 

NPL 

NSTOP 

JJ 

XMPli 

VEL 

DELTAT 

DELTAL 

FINAME 

FNAME2 

FK'S 

YINPUT 

INDEX 


INPUT  PROFILE 

FEELER  THRESHOLD  HEIGHTS 

OUTPUT  DISPLACEMENTS  (EQUIVALENCED  TO  DRV1) 

OUTPUT  VELOCITIES  (EQUIVALENCED  TO  DRV2) 

OUTPUT  ACCELERATIONS  IN  IN./SEC2  (EQUIVALENCED  TO  DDRV2) 
OUTPUT  ACCELERATIONS  IN  g's  (ACCISS(I)/386) 

MAXIMUM  ACCELERATION 

MINIMUM  ACCELERATION 

SUM  OF  THE  ACCELERATIONS  SQUARED 

FMS  ACCELERATION 

BCD  VARIABLE  IDENTIFICATION 

BCD  INPUT  FILE  IDENTIFICATION 

BCD  CHARACTER  TO  RING  TELETYPE  BELL 

CONTROLS  TELETYPE  PRINTOUT 

PRINTOUT  TIME  INTERVAL  ON  TELETYPE 

REAL  TIME 

CONTROLS  FILE  PAGING 
CONTROLS  STOPPING  OF  PROGRAM 
END  OF  INPUT  FILE  BNDICATOR 
TANK  SPEED  IN  MPH 
TANK  SPEED  IN  IPS 
REAL  TIME  INCREMENT 

QUANTA  OF  LENGTH  BETWEEN  INPUT  PROFILE  POINTS 
NAME  OF  INPUT  PROFILE  FILE 
NAME  OF  OUTPUT  FILE 

TEMPORARY  DEPENDENT  VARIABLES  IN  RUNGE  KUTTA 
PROFILE  INPUT  POINT  (FROM  FILE) 

INDEX  COUNTER  FOR  RK 
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NOTE:  Many  of  the  variables  stored  in  array's  are  assigned  individual  names 
as  well  as  subscripts  by  use  of  the  "EQUIVALENCE"  statement.  This  allows  for 
the  use  of  either  individual  namea  or  subscripts  in  the  program  whichever 
is  more  convenient. 

EXAMPLE: 

COMMON  VNAME1,  VNAME2,  VNAME3 
DIMENSION  VN*ME(3) 

EQUIVALENCE  (VNAME1 ,  VNAMO(l)) 

RESULT: 

VNAME1  -  VNAME(l) 

VNAME2  -  VNAME(2) 

VNAME3  -  VNAME(3) 

The  following  variables  are  used  in  this  manner: 

FORCH,  FORCD,  FORCN,  FORCK,  SPDEF,  PSPDF,  DISPL,  VELCTY,  AND  ACCISS. 
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COMPUTER  PROGRAM  FOR  SIMULATING  DYNAMIC  RESPONSE 


TRUCK 


DICTIONARY  of  PROGRAM  VARIABLES 


Program 


OF  M37  TRUCK,  AND 


I SLIB, DIFFER 
2SLIB.ALGEBR 
3*RPC 
4SNDM 
5STTY , 120 

100  COMMON  NSTEPS  .DELTA T  ,A  ,B  .HASS  .MASS  1  .MASS2.CP0S  l  .CPOS 2,C»SEGI  ,CNEQ2 
110  COMMON  C  1 1  ,C21  .INRT1A  .FORCWl  .F0RCW2.SPDEF1  .SPDF.F2.DSPDFI  .DSPDF2 
120  COMMON  FORCX  1  ,F0RCK2,VAR 1 .VAR2.VAR3. VARA, VARS, VAR6.VAR7. VARS 
130  COMMON  ZETA I ,Z£TA2,DZETA2,ETA 1 ,El A2.DETA2.NU 1 .NU2.DNU2 

140  COMMON  MU  I .MU2.DMU2 

130  COMMON  THRESH  (&0)  .GAMMA  (20) ,SEQDEF(20) ,Y C47) 

160  REAL  NU 1 .NU2.MU 1 .MU2.NASS .MASS I .MASS2.INRTI A 

163  IBELL=458732 

170  PRINT  1 

180  1  FORMAT ( I8X  ,33(  1H*) ./, 18X  , lHtf.SlX  ,  1H*,/, 

1906  I8X.33H*  MURPHY  S  M-37  BACK-UP  MODEL  *,/, 

2006  18X,1H*,31X,IH*,/,18X  ,33(  1H ,///) 

21 OC . DATA  INITIALIZATION . — . 

220  A:64.8 

230  B=47.2 

240  INRTIA=23446. 

230  MASS =8. 03 

260  MASS 1=. 94 

270  MASS  2=  .8 

280  CPOS  1:11. 8 

290  CNEG1=23.8 

300  CP0S2=I0.2 

310  CNEG2:44. 

320  THRESH <1>=7. 02 

330  THRESH (2):4. 3 

340  THRESH<3)=2.34 

330  THRESH (4):. 81 

360  THRESH (3):. 09 

370  THRESH(6):.09 

380  THRESH (7):. 81 

390  THRESH (8): 2. 34 

400  THRESH (9): 4. 3 

410  THRESH (10) :7. 02 

420  THRESH (1 1 ):7.02 

430  THRESH(I2):4.3 

440  THRESH (13):2, 34 

430  THRESH(14):.B1 

460  THRESH (13): .09 

470  THRESH (1 6):. 09 

480  THRESH(!7):.81 

490  THRESH (18)=2.34 

300  THRESH (19): 4. 3 

310  THRESH (20): 7, 02 

520  GAMMA ( 1 ) =41 1.73 

530  GAMMA (2): 506. 25 
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TRUCK  CONTINUED 


MO  GAMMA  <3>=587. 25 
550  GAMMA <4 >=€48. 

560  GAMMA <5> =668. 25 

570  GAMMA <6> =668. 25 

580  GAMMA (7)=648. 

590  GAMMA (8 >=58 7. 25 

600  GAMMA <9> =506. 25 

610  GAMMA (1 0)=41 1 .75 

620  GAMMA ( 1 1 )  =  4 1 1  •  25 

630  GAMMA (12) =506. 25 

640  GAMMA (13>=587. 25 

650  GAMMA (14)=648. 

660  GAMMA (15)=668. 25 

670  GAMMA (16)=668. 25 

680  GAMMA (17)=64B. 

690  GAMMA (18)=587.25 

700  GAMMA  (19  >  =  506. 25 

710  GAMMA (20)=41 1.25 

720  ZETA l=-5. 189 

730  ZETA 2=0. 

740  ETA  Is. 028 

750  ETA  2=0. 

760  NU  1  =  - 1 .  08 

770  NU2=0. 

780  MU  U-  1.245 

790  MU  2=0. 

800  NSP=0 

810  POINT  1  =  0. 

820  P0INT2=0. 

830  SDZ2SQ  =  0. 

840  T=0. 

850  NSTOPsO 

860  NPL=1 2 

870  JJ= 1 

880  DO  10  1=1,47 

890  10  Y (I )=0. 

900C . DATA  READ  IN . - 

910  PRINT,"************TYPE  IN  THE  DATA  FOR  THE  FOLLOWING 

920  PRINT, 

930  PRINT, "EXTERNAL  INPUT  DATA 

940  PRINT , 

950  PRINT,  NAME  OF  PROFILE  INPUT  FILE 

960  READ  2.FINAME 

970  2  FORMAT (A  6) 

980  PRINT, "DA TA  VARIABLES" 

990  PRINT, 

1000  PRINT, ‘TRUCK  VELOCITY  IN  INCHES  PER  SECOND  (REAL) 

1010  READ,VEL 

1020  PRINT,  PROGRAM  VARIABLES 

1030  PRINT, 
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TRUCK  CONTINUED 


1040  PR I  NT,  "NUMBER  OF  SPACES  BETWEEN  PROFILE  POINTS  (INTEGER) 

1090  READ  ,NSPACE 

1060  PRINT, "NUMBER  OF  STEPS  IN  RKG  (INTEGER) 

1070  READ  ,NSTEP5 

1080  PRINT,  EXTERNAL  OUTPUT  DATA 

1090  PRINT, 

1100  PRINT,  NAME  OF  OUTPUT  DATA  FILE 

1110  READ  2.FNAME2 

1120  PRINT, 

1 1 30  PRINT, “********END  OF  DATA  INPUT *********“ 

1140  PRINT  3 

1130  XMPH:VEL/17.6 

1160  3  FORMAT (/////, "T  ",  l OX  ,*Y(  1 )"  ,7X  , "AZETA2" ,5X , 

1170*  "ARMS",/) 

1180  WRITE(2|6) 

1190  6  FORMAT<//,37X,45(lH*),/,37X,lH*,43X,IH*J/,37X, 

1200*  45H*  MURPHY  'S  M-37  TRUCK  PROGRAM  OUTPUT  FILE  *,/, 

1210*  37X , 1H*,43X ,1H*,/,37X,43(IH*) ,///) 

1220  WRITE<2|U)XMPH,NSTEPS 

1230  WRITE <2|9) 

1240  II  FORMAT (/,9HVEL0CITY  = ,F6.2, 15H  MILES  PER  HOUR.50X, 

1230*  3 5H NUMBER  OF  STEPS  IN  RKG  INTEGRATION: , I  4) 

1260  9  F0RHAT(//,13X,9(1H-),1 2HDISPLACENENT ,8<IH-),3X,10(1H-), 

1270*  8HVEL0CITY , 1 1 < 1H-)  ,3X,9C 1H-) , 1 2H ACCELERATION  , 

1280*  8( 1H-) ,2X ,7HC-G  RMS  ,/,X .4HTINE ,2X ,4HY ( I ) ,X , 

1290*  3<2X,3HC-G,4X,3HPITCH,3X,3HFR-AX,3X,5HRE-AX,2X), 

1300*  X , 3HACCEL ,//) 

1310  DELTAT :3,07/VEL 
1320  CALL  OPENFd  , FI  NAME) 

1330  4  FORMAT (E20> 10) 

1340  30  POINT  1=P0INT2 

1350  IF(JJ-2)55, 40.33 

1360  33  READ(1,4)P0INT2 

1370  CALL  EOFTST (  1  ,JJ ) 

1380  GO  T0(60,40)JJ 
1390  40  P0INT2=P0INT1 

1400  NST0P:NST0P+1 
1410  NSPACE: 1 

1420  60  PSTEP:(POINT 2-POINT  1 ) /NSPACE 

1430  NSP=0 

1440  100  1:46 

1450  70  Y (I +1 )=Y (I ) 

1460  1=1-1 

1470  IFCI )70,75,70 

1480  75  Y ( 1 ):POINT  1 

1490  POINT l:POINT l+PSTEP 

1500  NSP:NSP+1 

1510  CALL  DIFFEQ 

1520  SDZ  2SQ:SDZ  2SQ+DZETA 2*DZETA2 

1530  RNSDZ 2:SQRT ( (SDZ  2SQ *DELTAT )/T ) 
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TRUCK  CONTINUED 


1540  5  FORMAT (6 (X, 610. 4)) 

1550  NPL=NPL+I 

1560  AZETA2sDZETA 2/386. 

1570  ANU2:DNU 2/386. 

1580  AMU 2sDMU 2/386. 

1590  RMSAZ  2:RNSDZ 2/386* 

1600  VRITE(2|7)T  .POINT  I  ,ZETA  I  ,ETA I  ,MU  I  ,NU  I  .ZETA2.ETA2.MU 2.NU 2, 
161 OA  AZETA  2.DETA  2. AMU  2.ANU  2.RMSAZ2 

1620  7  FORMAT (F6.I.F5.1.1 3F 8.3) 

1630  IF(NPL-54)I20,I 10,120 

1640  110  WRITE(2|9) 

1650  PRINT  5.T  ,Y  <1 )  .AZETA2.RMSAZ2 
1660  NPL-0 

1670  120  IF(NSTOP-47)80,90,80 
1680  80  TsT+DELTAT 

1690  IFCNSP-NSPACE )1 00, 50, 100 
1700  90  CALL  CL0SEF(2,FNAME2,2) 

1702  PRINT  1 2, (I BELL ,1 =1,40) 

1704  12  FORMAT (40A I ) 

1710  PRINT  8.FNAME2 

1720  8  FORMATt///,  DETAILED  OUTPUT  IS  IN  FILE  ,X,A6,/> 

1750  CALL  EXIT 

1740  END 
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01 FFEQ 


100  SUBROUTINE  D1FFEQ 

110  REAL  NU 1 ,NU2,MU 1 ,MU2,MASS  .MASS  1  .MASS2.INRTIA 

120  H=DELTAT /NSTEPS 

130  RH-l./H 

140  I NDEX  =  0 

1 50  100  INDEX  =  INDEX+1 

160  VAR1=ZETA1 

170  VAR2=ZETA2 

180  VAR3=ETA  1 

190  VAR4=ETA2 

200  VAR5=NU1 

210  VAR6=NU2 

220  VAR  7rMU 1 

230  VAR8-MU2 

240  PZETA  2=Z£TA  2 

250  PETA 2= ETA  2 

260  PNU2=NU2 

270  PMU2=MU2 

280  CALL  ALGEBR 

290  FK11=H*VAR2 

300  FK  12=H*(1/MASS)*(F0RCK  1+C 1 1*DSPDF1+F0RCK 2+C21*DSPDF2-MASS*386.0) 

310  FK13=H*VAR4 

320  FK I 4=H  * ( 1 /I  NR  TI A ) * (A *FORCK 1+ A  *C  1 1 *DSP  OF  1 -B*FORCK 2-B  *C  2 1 *DSPDF  2) 

330  FK15=H*VAR6 

340  FK 16=H*( 1 /MASS  J )*(-FORCK 1-C 1 1 *DSPDF 1+FORCW l -MASS  1*386.0) 

350  FKI7=H*VAR8 

360  FK18=H*(l/MASS2)*(-FORCK2-C2l*DSPDF2+FORCV2-MASS2*386.0> 

370  VARUZETAI+FKIU.5 

380  VAR2=ZETA2+FKI2*.5 

390  VAR 3=ETA  1+FK 134.5 

400  VAR4=ETA2+FK 14*. 5 

410  VAP5=NU1+FK 15*. 5 

420  VAR6=NU2+FK 16*. 5 

430  VAR  7=MU 1+FK 1 7*. 5 

440  VAR8=MU2+FK 18*. 5 

450  CALL  ALGEBR 

460  FK21=H*VAR2 

470  FK22=K*( I /MASS )*(FORCK 1+C 1 1*DSPDF1+F0RCK  2+C21*DSPDF2-MASS*386.0> 

480  FK23=H*VAR4 

49 C  FK24=H*(I/INRTIA)*(A*F0RCK l+A*C 1 I*DSPDF1 -B*FORCK 2-B«C  21 *DSPDF2) 

500  FK25=H*VAR6 

510  FK26:H*( 1/MASS  1 )*(-FORCK 1-C 1 1*DSPDF1+F0RCV1 -MASS  1*386.0) 

520  FK27=H*VAR8 

530  FK28:H*( I /MASS  2>*(-F0RCK  2-C  21 *DSPDF2+F0RCtf 2-MASS  2*386.0) 

540  VAR IsZETA l  +  .29289322*FK  21  +  . 2071 0678*FK 1 1 

550  VAR2=ZETA2+.29289322*FK22+.20710678*FK  12 

560  VAR3=ETA 1+. 29289322*FK 23+. 2071 0678*FK 13 

570  VAR4=ETA  2+. 29289322*FK  24* • 207 1 0678*FK 14 

580  VAR5=NUI+.29289322*FK25+.20710678*FK 15 

590  VAR  6=NU  2+ . 29  289322*FK  26* .207106  78*FK 1 6 
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DIFFEQ  CONTINUED 


600  VAR  7=MU 1+ . 29289322*FK  27+. 2071 0678*FK 1 7 
610  VAR  gsMU  2+ . 29  289 3  22*FK 28+ . 20  7106 78*FK18 

620  CALL  ALGEBR 

630  FK  31  =H*VAR  2 

640  FK32=H*(1/MASS)*<F0RCK 1+C1 1*DSPDFI+FORCK2+C21*&SPDF2-MASS*3P6.0> 

650  FX  33sH*VAR  4 

660  FK34=H*(1/INRTIA)*<A*F0RCK  l+A*Cll*DSPDFl-B*F0RCK2-B*C2t*DSPDF2> 

670  FK35=H*VAR6 

680  FK36sH*(l/MASS 1  )*<-FORCK 1-C 11 *DSPDFl+FORCWl-MASS 1*386.0) 

690  FK37=H*VAR8 

700  FK  38=H*<  1  /NASS  2)*  (-FORCK  2-C  21  *DSPDF2+F0RCW2-MASS  2*386.0 ) 

710  VAR  IsZETA 1-. 707 106 78*FK 21+1 • 7071 0678*FK31 

720  VAR  2=ZETA  2- . 707 1 0678*FK  22+1.7071 0678*FK  32 

730  VAR  3:  ETA 1-.70710678*FK  23+1 . 7071 0678+FX 33 
740  VAR  4=ETA  2-. 707 1 0678*FK  24+1.7071 0678*FK  34 

750  VAR5sNU 1-. 7071 0678+FK 25+1 »70710678*FK35 

760  VAR  6=NU  2- . 7071 0678*FK  26+1.7071 0678*FK  36 

770  VAR  7=NU 1 - . 707 1 0678*FK  27+1 . 707 1 0678*FK  37 
780  VAR8=MU2-.7071 0678+FK 28+1 • 70  71 0678*FK  38 

790  CALL  ALGEBR 

800  FK41=H*VAR2 

810  FK42:H*( 1/MASS  )*(FORCK  1+C  1 l*DSPDFI+FORCK  2+C21*DSPDF2-MASS*386.0) 

820  FK45=H*VAR  4 

830  FK 44±H*(  1/INRTIA  )*(A*FORCX  1+A*C  1 1*DSPDFI-B*F0RCK2-B«C21*DSPDF2> 

840  FK45=H*VAR6 

850  FK  46=H  *  C 1 /MASS  1 ) *  < -FORC  X 1 -C 1 1 *DSPDF 1 +FORC W 1 -MASS  1*386.0) 

860  FK47=H*VAR8 

870  FK4S=H*<1  /MASS  2)*<-F0RCK  2-C  21 +DSPDF2+F0RCW2-MASS 2*386.0) 

880  ZETA  1=ZETA 1+FK  1 1*. 166667+.09  763107*FK  21+.56903559*FK31+FK41*. 1661 

890  ZETA2=ZETA2+FK 12*. 1 66667+. 09 763 107*FK 22+. 56903559*FK 32+FK 42*. II 

900  ETA  1:ETA 1+FK  13*.  166667+.09763107*FK23+.56903559*FK33+FK43*. 1661 

9 1 0  ETA2=ETA  2+FK  1 4*. 166667+. 09763 107*FK  24+. 569 03559 *FK34+FK 44*. t  «6« 

920  NU1=NU1+FK 15*. 1 66667+. 09 7631 07*FK25+.56903559*FK35+FK 45*. 166667 

930  NU2:NU2+FK 1 S*. 1 66667+. 09 763 1 07*FK  26+. 56903559+FK36+FK 46*. 1 66667 

940  MU 1:MU 1+FK  17*. 1 66667+, 09 763 1 07*FK  27+,56903559*FK37+FK  47*. 1 66667 

950  MU2:MU2+FK 18*. 166667+. 09763 1 07*FK 28+. 56903559*FK38+FK 48*. 166667 

960  DZETA2=(ZETA2-PZET*2)*RH 

970  DETA  2: (ETA 2-PETA2)*RH 

980  DNU2=(NU2-PNU2)*RH 

990  DMU  2= (MU  2-PMU  2 ) *RH 

1000  IF (INDEX-NSTEPS  >100,200,200 

1010  200  RETURN 
1020  END 
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ALGEBR 


100  '  SUBROUTINE  ALGEBR 

110  REAL  NU 1 ,NU2,MU 1 ,MU2,HASS  ,MASS  1 , MASS 2,1 NRT1 A 

120  DO  100  1=1 , 10 

130  SEGDEF (I  )=Y(I  )-VAR5-THRESH <1  ) 

140  IF (SEGDEF (I >  )50, 1 00 , 1 00 
130  30  SEGDEF (I )=0. 

160  100  CONTINUE 

170  DO  200  1=11,20 

180  SEGDEF (I >:Y (I +26>-VAR 7-THRESH  (I ) 

190  IF (SEGDEF (I ) ) 1 30,200,200 
200  130  SEGDEF (I  ):0. 

210  200  CONTINUE 
220  FORCW|:0. 

230  F0RCW2:0. 

240  DO  300  I:  1, 10 

230  F0RCW1:F0RCV1+GAMMA (1 >*SEGDEF(I > 

260  300  CONTINUE 
270  DO  400  1:11,20 

280  F0RCW2:F0RCW2+GAMMA (1  )*SEGDEF(1 > 

290  400  CONTINUE 

300  SPDEF 1 =VAR 5-VAR 1-A*SI  N (VAR3) 

310  SPDEF2:VAR7-VAR1+B*SIN(VAR3> 

320  DSPDF  l:VAR6-VAR2-A*VAR4*C0S (VAR 3) 

330  DSPDF2:VARB-VAR2+B*VAR4*C0S(VAR3> 

340  FORCK 1: 19,54*(SPDEF 1**3)- 192. 424SPDEF1 ♦SPDEF 1+913.55*SPD£F 1 

350  FORCK 2: 1 ,39*(SPDEF2**3 )- I . 24*SPDEF2*SPDEF2+307. 72*SPDEF2 

360  IF (DSPDF I )600, 50 0,500 

370  500  C  1 1:CP0S  1 

380  GO  TO  700 

390  600  C  1 1  =C N EG  I 

400  700  IF (DSPDF2)800,900,900 

410  900  C2|:CP0S2 

420  RETURN 

430  800  C2I :CNEG2 

440  RETURN 

450  END 
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Variable 
NSTEPS 
D  ELTAT 
A 
B 

MASS 

MASSI 

MASS2 

CPOSI 

CP0S2 

CNEGI 

CNEG2 

N  SPACE 

VEL 

ZETAI 

ZETA2 

DZETA2 

AZETA2 

ETAI 

ETA2 

DETA2 

NUI 

NU2 

DNU2 

ANU2 

MUI 

MU2 

DMU2 

AMU2 

v 


Dictionary  of  Variables 

_ _ _ Des  cription _ 

NUMBER  OF  STEPS  IN  RUNGE  JCUTTA  SOLUTION 

REAL  TIME  INCREMENT 

DISTANCE  FROM  FRONT  AXLE  TO  CG 

DISTANCE  FROM  REAR  AXLE  TO  CG 

MASS  OF  VEHICLE 

MASS  OF  FRONT  AXLE 

MASS  OF  REAR  AXLE 

DAMPING  COEFFICIENT  FOR  POSITIVE  MOTION  OF  FRONT  AXLE 
DAMPING  COEFFICIENT  FOR  POSITIVE  MOTION  OF  REAR  AXLE 
DAMPING  COEFFICIENT  FOR  NEGATIVE  MOTION  OF  FRONT  AXLE 
DAMPING  COEFFICIENT  FOR  NEGATIVE  MOTION  OF  REAR  AXLE 
NO.  OF  INTERPOLATION  POINTS  BETWEEN  PROFILE  POINTS 
VELOCITY 

VERTICAL  DISPLACEMENT  OF  CG 
VERTICAL  VELOCITY  OF  CG 
VERTICAL  ACCELERATION  OF  CG 
VERTICAL  ACCELERATION  IN  G’S  OF  CG 
PITCH  DISPLACEMENT  ABOUT  CG 
PITCH  VELOCITY  ABOUT  CG 
PITCH  ACCELERATION  ABOUT  CG 
VERTICAL  DISPLACEMENT  OF  FRONT  AXLE 
VERTICAL  VELOCITY  OF  FRONT  AXLE 
VERTICAL  ACCELERATION  OF  FRONT  AXLE 
VERTICAL  ACCELERATION  IN  G’S  OF  FRONT  AXLE 
VERTICAL  DISPLACEMENT  OF  REAR  AXLE 
VERTICAL  VELOCITY  OF  REAR  AXLE 
VERTICAL  ACCELERATION  OF  REAR  AXLE 
VERTICAL  ACCELERATION  IN  G’S  OF  REAR  AXLE 
PROFILE  SHIFT  REGISTER 
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Variable 


Description 


POINT I 

POINT 2 

PSTEP 

SDZ2SQ 

RMSDZ2 

FINAME 

FNAME 

NSP 

JJ 

NSTOP 

I 

#  PROGRAM 


PAST  PROFILE  POINT 
PRESENT  PROFILE  POINT 

PROFILE  INTERPOLATION  INCREMENT  =[(P0INT2  -  POINTl)/NSPACE] 

SUM  OF  THE  SQUARES  OF  DZETA2 

RMS  OF  DZETA2  [RMS  (in. /sec2)]  RMSAZ2  RMS  (g)  OF  CG 
NAME  OF  PROFILE  FILE 
NAME  OF  OUTFUT  FILE 

PRESENT  STEP  IN  PROFILE  INTERPOLATION 
END  OF  PROFILE  FILE  INDICATOR 

STEP  COUNTER  FOR  T  AFTER  LAST  POINT  IS  READ  FROM  PROFILE  FILE* 
NO  SPECIAL  SIGNIFICANCE  USED  FOR  ALL  "DDloop"  COUNTS 

STOPS  AFTER  LAST  PROFILE  POINT  IS  READ  AND  SHIFTED  THROUGH  Y. 
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